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STAYLAM: A FORTRM PROGRAM FOR THE SUCTION 


TRANSITION ANALYSIS OF A YAWED WING 
LAMINAR BOUNDARY LAYER 

By James E. Carter 
SUMMARY 

A computer program called STAYLAM is presented for the computation of the 
compressihle laminar houndary-layer flov over a yawed infinite wing including 
distrihuted suction. This program is restricted to the transonic speed range 
or less due to the approximate treatment of the compress ihility effects. The 
prescrihed suction distribution is permitted to change dis continuously along 
the chord measured perpendiciolar to the wing leading-edge. Estimates of 
transition are made by considering leading-edge contamination, cross-flow 
instability, and instability of the Tollmien-Schlichting type. A program 
listing is given in addition to user instructions and a sample case. 

- ; : INTRODUCTION 

At the present time there is significant effort being made to implement 
boundary-layer suction on a wing to maintain laminar flow thereby resulting 
in a net drag reduction. Clearly such studies req[uire a computer program 
which analyses the compressible laminar boxindary layer on a swept wing and 
includes tests based on the latest technology to determine whether or not 
transition will occur for a given suction distribution. Based on the current 
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state of the art for houndary-layer computations for finite swept wings and 
that for transition estimates, it is clear that such a program would he most 
complex and thus difficult to use. Hence, in the present program a nxmher of 
approximations have been made in order to simplify the analysis; nonetheless, 
this program should be useful, particularly for preliminary design. 

SYMBOLS 



P 

Q 

Re 

T 


u 


U 

V 


w 

X 


speed of sound 
pressure coefficient 
skin- friction coefficient 
reference length 
Mach number 

indices for x- and z-directions , respectively 

static pressure 

free stream velocity 

free stream Reynolds number 

static temperature 

velocity component in x-direction 

velocity component in direction of inviscid streamline 
cross flow velocity .component: 

transformed velocity at the boundary-layer edge in the x-direction. 
velocity component at the boundary-layer edge in the y-direction 
(V = V' = sin i(j) . ■ 

velocity component in the z-direction 

coordinate along the surface measured perpendicularly to the leading 


edge 
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y 

z 

a 

Y 
Ax 
Az 

6 * 

0 

^2 

y 

V 
p 

T 

(}> 


coordinate along the surface measured parallel to the leading edge 

coordinate measured perpendicularly to surface 

weighting factor in finite-difference scheme 

ratio of specific heats 

grid spacing in x-direction 

grid spacing in z-direction 

displacement thickness 

momentum thickness 

q 2 du* 

pressure gradient parameter = — 
molecular viscosity coefficient 
kinematic viscosity coefficient 
density 

shear stress at surface 

angle hetween direction of flow at the houndary-layer edge and the 
X-direction 
angle of shear of wing 


Superscripts : 

* dimensional, untransformed q.uantity 

' dimensional quantity after Stewartson compressihility transformation 


Subscripts: 

e edge of boundary layer 

n normal to leading edge Of wing 

X in the x-direction 

y in the y-direction 

0° free stream quantity 
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GENERAL DESCRIPTION 


The present program, STAYLAM, was developed by modifying a program 
presented by Beasley (ref. l) for the calGiilation of the incompressible laminar 
boundary layer and prediction of transition on an infinite sheared wing. In 
Beasley ' s program the second-order accurate Crank-Nicolson finite-difference 
scheme is used to compute the boundary layer from the attachment line to some 
desired point downstream. These boundary-layer results are then analyzed to 
determine whether or not leading-edge instability or cross-flow instability 
occurs. The Owen-Randall criterion is used for the cross-flow instability 
test. The Tollmien-Schlichting type of instability is estimated by using a 
correlation given by Stuart (ref. 2) of the critical Reynolds number as a 
fiuiction of the external pressiire gradient. The point where transition is 
completed is then estimated by using a correlation given by Granville (ref, 3). 
These same tests are used in the present program; the only modification which 
has been made is that the input quantities to these tests are the actual 
compressible values, not the corresponding incompressible values given by the 
Stewart son transformation. It should be noted that the transition tests in 
the present program can be replaced or supplemented with relative ease. 

The Beasley program has been modified by the inclusion of distributed 
wall suction, compressibility effects, and the finite-difference scheme has 
been generalized to be of arbitrary accuracy between first and second order 
in the streamwise, marching variable. Figure 1 gives a typical distribution 
of suction velocities and explains the nomenclature used in inputting this 
distribution to the program. Note that the suction is allowed to change 
dis continuously at a prescribed nimiber of locations along the airfoil. It 
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was found that the Beasley program, modified to include suction, gave distri- 
butions of the local skin-friction coefficient which showed significant 
oscillations in a region of discontinuous suction. These oscillations were 
eliminated by using a first-order accurate finite-difference scheme in the 
streamwise marching variable, instead of the second-order accurate Crank- 
Hicolson scheme used by Beasley. The oscillations were expected in using 
the Crank-Bicolson scheme due to its known neutral stability in the wall 
region. The first-order accurate scheme suppresses the oscillations caused by 
the discontinuous suction since it has greater damping. For the same accuracy 
the first-order scheme requires more streamwise grid points than the 
second-order scheme; however, calcijlations showed that both schemes yield about 
the same result if approximately 100 grid points are used from the leading to 
trailing edge. Appendix A gives further details of the finite-difference 
s cheme . 

The Stewartson (ref. 1;) transformation has been used to account for 
compressibility effects. Details of this transfoarmat ion along with the 
Blasius transformation used by Beasley and the coordinate system are presented 
in Appendix B. The Prandtl number is assumed to be unity and the total 
temperature is assumed to be the same as the free stream value. A further 
approximation is made in the treatment of the streamwise pressure gradient 
term which allows the incompressible infinite sw‘ept wing equations to be 
obtained after the Stewartson transformation. This latter approximation 
restricts the use of the present program to speeds in the transonic range or 
lower. Since the present interest in LFC (laminar flow control) is in the 
transonic speed range it is felt that this simplifying approximation is 
justified. 




Several calculations were made to verify the accuracy of STAYLAM. First, 
the incompress ihle boundary layer on a circular cylinder with a constant suc- 
tion velocity was computed and comparisons were made with the results obtained 
by Terrill (ref. 5) for the same case. Excellent agreement was found in the 
momentum and displacement thicknesses, and skin friction distributions. The 
estimate of the separation point from the present program differed from that 
of Terrill by only 0.1°. A second test case was a comparison between the 
results from the present program and the analytic solution for the 
compressible asymptotic suction profile (ref. 3) which is obtained by applying 
constant suction on a flat plate. Excellent agreement was obtained in the 
streamwise and normal velocity distributions. 

PROGRAM USAGE 

The program was written in the FORTRAN programming language for use on 
the Contro]. Data 6000 Series Computer Systems under the NOS. 1.1 operating 
system at Langley Research Center. Included in the output is a plot 
of the distributions of the x and y skin-friction coefficients, 0^ and 
C^ , versus the non-dimensional surface distance measiored perpendiciolar to 
the leading edge. Some modification to the program might be required to 
obtain plots on a different computer system. 

The input and output for STAYLAM are discussed in the next two sections. 
The program listing is given in Appendix D. and a sample case is discussed in 
Appendix E. 



INPUT 


Read 

Order 

Variables 

Format 

1 

ITS, J, TOL, m, DS, USTEP, ^ 

2I5,5F10.5 

2 

IBLC 


15 

3 

MSMAX ' j 

(Skip if 

15 

h 

wall(ms), sds(ms) 

^ IBLC = 0) 



MS = 1, MSMAX 


2F1.0.5 

5 

INC, AMINE3D, GAMMA 


I5,2F10.2 

6 

B 


8aio 

T 

INT3, GH 


15, Flo. 5 

8 

ISY, INT4, RHO '1 (Skip if INT3 = O) 

,2I5,F10.5 

9 

XA(N), ZA(n) j 


2E16.8 


(INT4 pairs of valr.es) 


10 

ISP 


15 

11 

IFPT, nrai, NLIST, DX 


3I5,2F10.5 

12 

OPX(N), N= 2, NLlSTPl ( 

;Skip if IFPT ^ 2) 

8F10.5 

13 

c 


8A10 

li^ 

PSI, DTRIP 


8F10.5 

15 

INTRL, IFR 


215 

l6 

RNL(i), I = 1, INTRL 


8F10.5 

IT 

INTV, L 


215 

18 

XV(I), UM(I), I = 2, LPl 

2E16.8 

19 

MGRAD , XV(:i, \ . (if INC = 

= 0) 

2F10.5 

20 

VGRAD, XV(1) (if INC = 

1) 

2F10.5 


Input which has been added to original Beasley program is underlined. 



21 

22 

The 

ITS 

J 

TOL 

DZ 

DS 

USTEP 

WP 

IBLC 

MSl/IAX 

WWALL 

SDS 


INC 

8 


ALPHA, SI, S3 (skip if MGRAD or VGEAL > O; 
KUE 


3F10.5 


definitions of these input variables are as follows ; 

Maximim number of iterations in subroutine ¥UV to calculate the u 
profile. 

Humber of steps in z direction. 

Iterative tolerance. 

Step length in z direction. 

Standard step length in x direction, made non-dimensional with 
’ respect to airfoil chord. 

Maximum permissible' increase in velocity at edge of boimdary layer 
across one step. 

Accuracy control on finite-difference expression for x derivatives 
WE = 0.5 gives second-order accuracy (Crank-Hicolson scheme); 

WF = 1. 0 gives first-order accuracy. 

Suction parameter. IBLC = 0, no suction; IBLC = 1, distributed 
suction permitted. 

Total n-umber of values of suction velo city. 

- 


Array of suction velocity values = ^ jj- . For suction 

^00 ' 

WWALL(MS) <0. 


Array of nondimensional locations along airfoil measured from 
attachment line, where value of suction velocity changes 
dis continuously,. Set SDS(MSIiAX) > surface distance from 
attachment line to trailing edge. 

Compressibility parameter. IHC = 1, flow i$ incompressible. 
Iirc = 0, flow is compressible and Stewartson transformation 



is used. 


AMINF3D 

GAMMA 

B 

INT3 


CH 

ISY 

iwt4 


EHO 


XA 


ZA 


Free Stream Mach numher. 

Ratio of specific heats; usually Y == 1-k. 

Main title . 

= 0: velocity data vill be given at x co-ordinates, that is, 

at distances from the attachment line measured aromd the 
airfoil surface. 

= 1: velocity data will be given at chord-wise stations. 

Airfoil chord, measured perpendicularly to leading edge.. 

= 0: airfoil is cambered. 

= 1: airfoil is symmetrical. 

Initially is the number of pa3.rs of co-ordinates to be read, 
subsequently becomes the total number of pairs of airfoil 
co-ordinates stored, including the lower surface and leading edge 

Rose radius of the airfoil in plane perpendicular to the leading 
edge. . 

Coordinates of geometric data, measured perpendicularly to the 
leading edge and in the plane through the leading and trailing 
edges . 

Coordinates of geometric data, measured perpendicularly to the 
plane through the leading and trailing edges . For a symmetrical 
airfoil XA and ZA are read from the leading edge to the trailing 
edge and include values at both points. For a caiibered airfoil 
the geometrical data are read from the trailing edge on the lower 
surface to the trailing edge on the upper surface. Values of XA 
for the lower surface must have negative signs, and values of 
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ISP 

IFPT 


INTI 

NLIST 

DX 

OPX 

C 

PSI 

DTRIP 

INTEL 

IPR 


ZA must have signs as appropriate. 

Surface parameter indicator. ISP = 2, upper surface calculation; 

ISP = 0, lower surface calculation. 

= 1: Complete print-out at end of all steps. 

= 2: Complete print-out at points given in list. 

= 3: Complete print-out at points at which Velocity data is given. 

= k: Complete print-out at points DX apart, where DX is given as 

data. 


In current program set IFPT = 3 and the print-out has heen 

modified so that the complete print-out (houndary-layer prof iles ) 
is printed every 10^ chord. This modification can he eliminated 
by several program changes in subroutine PRINT. 

Velocity profiles are printed out at values of z corresponding to 
n = 1,2,3,. . .(INTI + 1), 2(INT1) + 1, 3(lNTl) + 1,...,F 

Number of points in output list. 

Interval between listed output points (when IFPT = A), made 
non-dimensional with respect to the airfoil chord. 

Points at which full output is required. 

Sub-title. 

Angle of shear. 

Trip-wire diameter. 

Number of values of Reynolds number to be read. 

Q^L*: sec 

OO I 

= 1. Data Reynolds number = — ' 

00 

= 2, Data Reynolds number = — 


= 3. Data Reynolds number - 


cos it 
V«- 
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In present protraia use IFR = 2 as this Reynolds number definition 


RNL 

IKTV 


L 

XV 

UM 


MGRAD 


has been assumed in the suction velocity and in the skin-friction 
coefficient calculations. 

Reynolds number . 

= 1: Velocity data is given as U. 

= 2: Velocity data is given as U sec Tp. 

=3: Velocity data is given as C . 

IP 

If flow is compressible set INTV = 1. See UM description for 
further explanation of input in compressible case. 

Initially is the number of velocity data points read in, 

subsequently is the total number of points at which the velocity 
distribution is defined, including the attachment line. 

Coordinates of velocity data. Use the same sign convention as 
that for XA to indicate whether the upper or lower surface is 
to be computed. 

For incompressible flow, UM is initially the velocity data, 
subsequently, U. For compressible flow UM is initially M 

n 

After the Stewartson transformation, it is U. The velocity 
(or Mach number) data, XV and UM, are read from the attachment 
line towards the trailing edge, but attachment-line values 
must not be included and the data need not extend all the way 
to the trailing edge. 

Mach number gradient (nondimensioual) at the attachment line in a 

dM„ 

plane perpendicular to the leading edge = 

11 


VG-RAD Velocit;f gradient (nondimensional) at the attachment line in a 

<au'/Q» 

plane perpendicular to the leading edge = — sec ^ 
XV(l) Location of attachment line. Use same sign convention as that 

used for XA. 


Incidence of airfoil in streamwise plane. 

Quantities used in equations (53-5H) in reference 1. 

Parameter for more data on step. 

= 1: Read more data from ITS, 

= 2: Read more data from B. 

= 3: Read more data from ISY. 

= Read more data from IFPT. 

=5- Read more data from C. 

= 6: Read more data from IRTV. 

= J: Stop. 

The input is printed and laheled as described above. In addition the 
following quantities are also printed along with the input. 

STH Kon-dimensional surface distance measured from lower surface 

trailing edge. 

TH Transformed chordwise station X. Lower surfaces 0 = cos ^(2[x| 

upper surface, 6 = 2it - cos ^(2X - l). 

FSTH Second derivative of siu'face distance with respect to 6. Used 


ALPHA 
SI, S3 
KUE 


in cubic spline interpolation. 

FZTH Second derivation of Z (measured perpendicularly to the plane 

through the leading and trailing edges) with respect to 9. 
Used in cubic spline interpolation. 



sxv 


SXVINC 

U 

THXV 

FUTH 

FSVSINC 


Non-dimensional surface distance measured from attachment line 

to point at which the inviscid velocity, Mach number, or pressure 
coefficient data is prescribed. 

Transformed surface distance corresponding to SXV; same as x'/L* 

in equation (bU). „ 

a* u* e 
“ e n 

Non-dimensional transformed surface velocity = — 

Q ^OO OO 

Value of 6 at the points where the velocity data is prescribed, 
XV. 

Second derivative of U with respect to 6. Used in cubic 
spline interpolation. 

Second derivative of SXV with respect to SXVING. Used in cubic 
spline interpolation. 


OUTPUT 


The displacement and momentum thicknesses, skin-friction values, and 
results of the transition estimates are printed at each location on the 
airfoil where a boundary-layer calculation is made. These locations are 
determined by the USTEP criterion, or if this is satisfied, then the 
computation is made at regular DS intervals. Furthermore, since IFPT = 3, 
computations are also made at the same locations at which the velocity (or 
Mach number) data is prescribed. Note that these latter computations are 
only temporary; hence results from these stations do not form upstream 
conditions for the next dovmstream station. The skin-friction coefficients 
are not computed at these intermediate stations. 
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The hoiindary-layer profiles are printed at approximately every 10 per 
cent chord. The information printed at each houndary-layer station and an 


explanation of 

X 

S 

SCOMP 

U 

AME3D 

DU/D(S/L) 

DELTAl 

THETAl 

(DU/DZ)Z=0 

(DV/DZ)Z=0 
AIEFOIL SLOPE 


the "boundary-layer profiles is given as follows: 
Non-dimensional chord location. 

Transformed incompressible coordinate measured along 
airfoil from attachment line; same as x in equation 
(12) in Appendix B. 

Non-dimensional surface distance measured from attachment 
line. 

Non-dimensional transformed velocity at the boundary- 


layer edge = 


a^' u* 

00 Q 


M. 


-n 


a* Q* M 

Q ^co a 


Mach number at the boundary-layer edge. See equation (18) 
in Appendix B. 


Transformed invascid velocity gradient , 


dx * 




Scaled displacement thiclmess = . 

/uRe ' Q 

Scaled momentum thickness = / ■ \ . 

Scaled skin-friction coefficient in x direction = 
2y-l 

f Y-1 2 \ 

1 + ^ It \ o 

2 ) 9u | 

3z' 


1 + 

2 e 


z=0 


Scaled skin- friction coefficient in y direction = 
3 y -1 


1 + 


2(y-1) 


2 j. Xli. 

2 “e 


I 


3v | 

3zl 


z-0 


Local airfoil slope in degrees in plane perpendicxilar 
to the leading edge. 
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(DIMENSIONAL 

CHORD 

DELTAl/C 

THETAl/C 

CFX 

CFY 

CDFX 

CDFXINF 

CHlCOWEN- 

RANDALL) 

RTHETA 

RTHETCRIT 

LAI42 


Z)/ Multiplicative factor to convert scaled displacement 

and momentum thicknesses into actual thickness divided 
hy the chord. 

6«/L« 

e/L^^ 

Skin-friction coefficient x direction = 


X 


2U- 


5/2 


1 2 
2 ^oo^co 


DUDZ 


i/xRe 


z=0’ 


Skin-friction coefficient in y direction = 


- — p = 2 sin i|j 


■^DVDZl 
xRe z=0 


Skin-friction drag coefficient i:a direction perpendicular 
to leading edge hased on chord nieasured pex’peridicular 
to leading edge. 

Total skin- friction drag coefficient in free stream 
direction hased on chord measured parallel to free 
stream. 


Cross-flow Reynolds numher = X = Re. 


J ^oo 


d( z*/L* ) 


u*0 

e 


Reynolds nuraher hased on momentum thickness = — 

min 

where V*. is the minimum of V* or V*. 

min oj e 

..Critical momentum thickness Reynolds numher from Stuart 
(ref, 2) for prediction of instability in the 
Tolirniien-Schlichting sense, 

'2 du* 

_ , 0 s 

Pressure gradient parameter , ~ — ' d^ 

min 
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INSTAB. 

RTC-RTI 


LAM2BAE 

Z 

ZCOMP 

W 

U 

V 

STV 

CFV 

T 

RHOD 

l6 


= (THETAl)' 


2 00 j ' 2 °° ^ 


X dU 


a + 


(l,I|l^u2)3/2 Udx 


RE . NO. Momentiom thickness Reynolds niamber at estimated point of 

iCTinar instability as determined from Stuart correlation. 

Critical momentum thickness Reynolds number minus momentxim 
thickness Reynolds number at the point of laminar 
instability from Granville correlation. 


*^xi 


Average value of Xg = 


""tr /l + ^ '\ffeT 
A 2 “ 


X* 

tr 


X* 

1 

L* 


Transformed coordinate normal to surfac e, z . 

me 

Scaled coordinate normal to surface = - 


00 


X L* 

Transformed velocity component normal to surface, w. See 


equat i on ( 12 ) in Appendix B . 


u’ 


Velocity tangent to surface in x direction = 


u 

u • 


Velocity tangent to surface in y direction - ~ • 

^e Q* sin il; 

Velocity component in direction of inviscid 

^ U« 


streamline = 


u* 

s 


See Appendix C for further explanation 


Velocity component perpendicillar to direction of inviscid 

u* 

streamline (cross-flow component) = • See Appendix 


C for further explanation. 


T* ' 


Static temperature ratio = ^ 
Density ratio - . 

Prn 



APPENDIX A 


FINITE-DIFFEEENCE SCHEME 

The computational molecule for the finite-difference scheme is shown in 
the accompanying sketch. The point of evaluation moves from the midpoint 
hetveen lines m and m + 1 as a, the weighting factor, varies from 0.5 to 
1.0. With derivatives evaluated at the point x = (m + a)Ax, z = nAz the 
following finite-difference approximations are given. 

9u ^m+l,n ^m,n , , 9^u . 2% / \ 

_ = ___i_ _ji_ + _ i/2)Ax — ^ + 0(Ax ) (AT) 

3x 


3u v/’‘^m,n+l ^m,n-l\ . Y^m+l,n+l ^m+l,n-l^ . 2> 

^ = (1 - a)(— ‘—y * ^ 2Az ' > ' 


(A2) 


u 


u , T - 2u + u T u ^ 

:= (1 _ m,n-M m,n m,n-l ^ ^ m+l,n+l 


3z‘^ 


Az 


izu u , ^ T 

m+l,n m+l,n-l ', 

2 ^ 


Az 


(A3) 


-!* O(Az^) 


From, the truncation error it is seen that if a = 0.5 the scheme is second- 
order accurate, which is the Crank-Nicolson scheme used hy Beasley. First- 
order accuracy is obtained if a = T and this scheme is sometimes referred 
to as a fully implicit finite-difference scheme. In the program the 
weighting factor a is designated as WF:, and should be restricted to the 
range 0.5 < WF <1,0. It should be noted that if a = 0.5 the normal 


component of velocity which is printed is the value at m + 1/2, n; 
otherwise, for, a > 0 . 5 the value is at m + l,n. 





APPENDIX B 


GOVEBNING EQUATIONS 

The compressible, laminar boundary-layer equations are given as follows 
for the flow over an infinite, swept wing: 



Figure B-1. Yawed wing coordinate system. 
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The Stewartson transformation, which is given hy 


■p* a^ 
P* 

•^00 00 


z* a"' „ 

dz^ 

a'^ p* 

00 '^00 


^00^00 / V 3z' «’ 

a«p« '' 3x« p« ' 


is applied to equations (b 1) - (B3). In addition it is assumed that the 
Prandtl nimiher is unity, the total temperature equals the free stream value, 
and the viscosity coefficient varies linearly with the temperature. The 
transformed equations are: 


9u' ^ 8w' - 

aP"" a?" ° 


Y -1 t ,*2 . 2 , 

.811' , V 

1 + cos^ij 


. 2 , / e dx’ 
Q* sxn ?jj/ 


(BIO), 



The coefficient of the external velocity gradient in equation (b 9) gives rise 
to a coupling "between the x and y momentum equations which is not present in 
incompressible flow. This coefficient varies from its maximum value at the 
surface to unity at the botindary-layer edge as v* approaches its edge 
value Q* sin tJ;. This maximum value increases as the Mach number and sweep 
angle increase; nevertheless, this coefficient remains close to unity for 
flows in the transonic speed range, which is the present area of interest. 

For example the maximum value of this coefficient is 1.06 or less for ijj = 35° 
and free stream Mach numbers of one or less. In the present program this 
coefficient is set equal to unity and the x momentum equation becomes 







(Bll) 


Thus it is seen that with the given assumptions the Stewart son transformation 
converts the compressible equations into an equivalent incompressible 
formulation. This formulation is the starting point for Beasley's analysis 
which will be repeated for convenience. 

Beasley applied the Blasius transformation given by 



(B12) 


to equations (b 8), (Bll), and (BlO), respectively, and obtained 



X 
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9u ^ X dU ^ z,x dU - V 3u ^ 8v * 


9u , uz /X dU T \ T 3u 

t"* — ‘uS- ^ 


( 1 .^ 2 ) 2.®.* La 

^ ^ U dx ' 2 

3z 


(Bl4) 


3v ^ r ^ uz /X dU T ,,T 3v 

[v + — (jj— -1)]^ 


3^v 

3z" 


(B15) 


The surface 'bOTindary conditions including suction are given as 


z = Oj u=v= 0 




w = 




Q* 

^CO 


(Bl6) 


The edge conditions are 


Z -> 00 


u 1 , 


V 1 


(BIT) 


Hote that in equation (bi 6) the Mach number of the inviscid flow is given as 


m2 

M = 
e 




n ^ T-1 t ,2 2 , 

1 4- -^7-— H cos Ip 

2 00 T 


(B18) 


>r 


*/Rs 


The scaled, non-dimensional suction velocity. 


Q* 

^00 


is referred to as 


VrWALL in STAYLAM and is part of the input. 

After the solution is obtained in terms of transformed, incompressi'ble 
variables the corresponding compressible quantities are obtained as follows: 


22 


X /i + \2 (y-1) 


xi r A ^ ^ '4 \ 
"*'•4 V^^nl) 


Y-1 2 \ 

1 + ^ M \ 2(y-1) 


1 + 

2 °° 


= 1 + 


< {‘ ' fef • (s9 



/i + 

u* 

I 

Q* 

^co 

u + ^ 


\ 2 

V* 

V sin \jj 

Q* " 


Y 

y -1 

1 + M \ T* 

2 CX> \ 00 


Vl + ^ 

\ 2 e 


The moment-um and displacement thicknesses are given hy 


L* U , 1=1 ,,2 

\ 2 “ 


1 + j 42 \ 2 (y~ 1 ) 


u(l - u)dz 


^ = Ai\ - ; 

.» \l*1 4 


i + 2 ^m 2 \a(Y-D 
2 e 1 


\ 2 °o 


where the subscript e denotes the edge of the boundary layer. The local 
skin-friction coefficients are 


1 



2y-1 

Y-1 


(B2T) 




(B28) 




APPENDIX C 


VELOCITY COMPONENTS 


The resolution of the velocity components at any point in the bomdary 
layer into the strearawise and crossflow components, respectively, is given 
as follows: 


= u* cos (|) + V* sin (j) 


u* = V* cos (j) - u* sin $ 


These components are non-dimensionalized by the free stream velocity, Q*, 
and the Stewart son transformation is incorporated to give 


u* / 1 + -h” M 
s 2 “ 




Q* \l + 

\ 2 e 


u U cos ^ + V sin 4* sin 


/l + 

—H. = u[ 

Q* 

\ 2 e 


(v - u) sin cj) 


Note that at the boundary-layer edge 


A + Zzi \ 

.21 2 \ . . 2, 


u* = 0 
c 



APPENDIX D 


PROGRAM LISTING 

PPOGRAf^ STAyLANntgpUT»DUTPUT»TAPEl»If^PUT»T*PE2»OUTPUT) 1 

2 

r master laminar rcujndary layer 3 

COMMOM/SUPWUV/J, ITS, THL, T ECOM, 05,(56 1 # nZ,WF,SH,UMGH,ACi,CGN3,C0N6* 5 

1 rnN7,G, ANCnnI ,BD0C170) ,BNC170) ,BNO( 1705 ,BNDn70),CNn70) ,DDOC170) 6 

2,n>tn70 5 ,RD1,DN0{ 5 70 3 , D^D f 1 70 ) , IJM J ( 1 70 ) , UMG ( 1 70 3 , WM ( 1 70 5 , WH i ( 1 70 ) , 7 

3*-'Nn( 1 70 3 ,ANn( 1 703 6 

CCMMON/SFX/STWCi*,?) ,TH{5fc5),F8TH(3B53,TN7fl,EZTH(3(!>B) R 

C(3MMOM/ySANDU/UMC3A53 ,THXV(5653 . FUTHOfeS 3 , XV C565 3 , CPUM (3<»5) , 10 

15XV C565 3 ,SXVINC (36S) ,FSVS1NC(36B3 , L , 8 A TT , IMT3. CH , ISP 11 

r.O.MMON/RESULTS/WMPfUOD J 70 3 ,Vf^2( j 70) ,DELTA1 ,THETAl,NO,DUr)Z 12 

l,ovnz 13 

CfiMHON/REnM/XA C3653 *ZAf 3653 1« 

DIMENSION XAHLD(3653 ,7AHLDC365) 15 

COMHnN/TEST/BNl.aO) , TNTRL, IFR 16 

r,0*^-H0W/0PLlST/DPX C200 3 ,OP8{2003 ’ 17 

CCMMONycONPRES/IMC, AM1NF3D, AHE5D,8INP,CnSP,6AMMA,GAMi,GAM2,GAM3, 16 

ir,AMii,AMFS30, RATIO, SCOMP,ZCOHPCl703 ,Tn7O),RHRr)C170) 19 

REAL, HGRAD 20 

filMENSlON S(2),B(8),C(63 21 

niMENSlON UMIHU 703 ,UH2HC1703 ,VM2HC1703 22 

dimension wt^ALL US) ,SDS( 153 - 23 

DIMENSION VCFX2t?00 3 ,VCFy2(200) ,HSCOMPt2003 ,LAaF-L(83 29 

1 ,LA8C83 , 1NFO(503 25 

REAL HMAX( 10) 26 

27 

PI«3.1 U 159265 26 

DTRb, 01 7953292 29 

call pseudo 30 

CALL LEROY 31 

32 

C HEAD BOUMDAPY layer CALCULATION PARAMETERS S3 

r maximum NO, OF iterations, no, of steps IN z direction, iterative 39 

t TOLEPANtrE, STEPLFNGTHS IN Z AND X DIRECTIONS, AND MAX VARIATION IN U 35 

C ACROSS ONE step, weighting FACTOR FOR FINITE OIFFERFUCE SCHEME 36 

37 

CALL JPARAMSaNFO) 38 

ENCODF(80,/i25,LABa) 3 INFO (1 ) , I NFO C25) , INFO f 22 ) 39 

925 FORMATCA7,3X,2A10, 50X3 90 

29 READtt,6)IT8,J,T0L,n2,DS,USTEP,WF 91 

6 FORMAT(Pl5p5F10.55 hZ 

93 

WRlTEC2,i033 99 

103 F0RMAT(lHt/lX65(2H*»3// 95 

199X,*INPU7 + SOME COMPUTED QUANTITIES* 96 

2//1X65C2H*.)/) 97 

WRITEC2, 900)118, J,TOL,DZ, OS, U8TEP,WF 98 

«nn format C5H IT$a, 15, 2X,2HJ«,I5,2X,«HTOL»,F10.5,2X,3 HDZb,F10.5,2X, 99 

13HP8»,F10,5,2X,6HU8TEPii,F10,5,2Xi3HtxF»Fl0,8) 50 

51 

C READ BOUNDARY LAYER SUCTION OR INJECTION PARAMETER, IBLC. 52 

C IBLC*0, NO SUCTION OR INJECTION PERMITTED. 53 

C IBLC«T> DISCONTINUOUS SUCTION OR INJECTION PERMITTED, 59 

55 

READd , 77 3 IBLC 56 

77 F0RMATCI5) 57 

IF(IBLC.EQ.0)OO TO 76 58 

59 

C READ SUCTION OR INJECTION VELOCITY P AR AMETER , WW ALL (M8)R W/Q*6QRT (Ot/NU 60 

■ 27 


05) POOR QUALITY 



noo 


C ) which is ALtOWED TO CHANGE 0 1 SCONTINUOUSL, Y AT SnS(M8} LOCATIONS, A 61 

C TOTAL OF MSHAX LOCATIONS A»E PERMITTED, 62 

63 

READC1,77)MSMAX 6« 

PEAOn ,502HHWALL(M81 ,8ns(M8) ,M8«1,MSMAX) 65 

VsRITE(2,70) 

19 EORMATC/a 1RLC»1 discontinuous suction or injection given below*) 67 
WRITE(2,fl03 MSMAX 68 

60 FORMAT C/lX29H8UCTTnN OR INJECTION VELOCTTY,5XIOHS LOCATION, 69 

l5X6HMSMAXaI53 . 70 

WRlTE(2,«n{WWALLCMSi,SOS(HS),M8»i,MSMAX) 71 

61 FORMATn(.)Xp9HWwALLC«S3,20X7HSnS*:M8)/C10XF10,5,l5XF10,5) ) 72 

GO TO 82 73 

7P fRITEtSjRa) 79 

3ii FOR-ATC/* IBLCnO NO SUCTION DR INJECTION PERMITTED*/) 75 

wisALLn)*O,0S 80S(n«l 00000, 76 

62 CONTINUE 77 

78 

READ COMPRESSIBILITY PARAMETER , INC , AND FRCESTREAM MACH N0.,AMINF3D 79 

INCal FLOW IS INCOMPRESSIBLE 80 

INC«0 FLOW IS compressible, COMPUTATION MADE IN STEHARTSON 81 

TRANSFORMED VARIABLES 82 

83 

READ 85, INC, AMINESn , gamma 89 

85 FOPMATCI5,2FT0.2) 85 

IFnNC,E0,l) AMTNF3DB0, 86 

PRINT B6»INC,AMINF3D,CaHNA 87 

86 FORMAT!/* INC»*T2,* AMlNF3D«*E16,ft,* GAMMA«*E1 6,81 88 

IF(INT.,EO,n r,n TO 86 69 

PRINT 87 PO 

87 FORHATf* compressible FLOWCINCpO), STEWARTSON TRANSFORMATION USED 91 

1*1 92 

GO TCI 90 93 

86 PRINT 89 99 

69 FORMATC* incompressible FLOW (INCbII,*! 95 

90 CONTINUE 96 

97 

CHMSTANTS needed for COMPRCSSieLE FLOW CALCULATION 98 

99 

Gm1»GAmma» 1, too 

GPlaGAMMA+l, lOl 

GAM)»i(5’,*CAMWA-l,l/t2.*r.Ml) 102 

GAMgBGPl/(g,*GMll 103 

SAMSbGAMMA/GNI 109 

GAM9«(2,*GAMMA-l,)/(GAMHA»l,) 105 

AMF83D*’,5*GM1*AMINF3D**2 106 

C read main TITLE (COLS l->721 107 

79 READ{1,1511B 108 

151 FOPMAT(8A103 109 

teRlTEC2,920ia 110 

920 FORMAT(1HO,3HP» ,8At03 111 

112 

C REAP geometry PARAMETER AMD CHORD LENGTH 113 

C IF INT3B0 VELOCITY data WILL RE GIVEN AT INTERVALS MEASURED AROUND THE ll9 

C AEROFOIL SURFACE AND NO GEOMETRICAL DATA IS REQUIRED 115 

116 

PEA0(1,72)1NT3,CH 117 

72 p0RmaT(15,F10,51 118 

WRITE(2,901 1INT3,CH 119 

901 FaRMAT{lH0,5HlNT3«,I5,2X,3HCHli,F10,53 120 
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121 

C ^’U8T PIST>Mr.E9 AROUMr SURFACE BE CALCULATED 122 

IE(It«T3,EQ,05GO TO 75 123 

12« 

C READ GEOK’ETRIC DATA 125 

C SYwHfTRY parameter* MO, OE PAIRS OF CO-ORDIMATE8* NOSE RADIUS. I2fc 

r IE ISYsO AEROFOIL IS CAMRfREO# IF 18Y«1 AFPOFOIL 18 SYMMETRICAL, 127 

128 

7fc PFADU , 1 00 3 ISY,INTa,RHO 129 

too format C2I5,2F1 0,5) 130 

MRITE(S,a02)TSY, IMT<!4,RMO JSl 

«02 FnRMATtlH0,UHTSy»,I5,?X,5HIMTA»,I5,2X,«HRH0B.E18,8) 132 

RHObRHO/CH 153 

13A 

C READ CIMTti) pairs OF X AND Z . 135 

139 

Nlsf IMT6-n*ISY + t 157 

M2B(lMTfl-l)*TSv+TMTU 138 

REAn(l,510HXAHLr'(:M),ZAHLD(M),N»Nl ,N21 139 

510 FORMAT (?Elh. 8) UO 

502 FDRMAT(2F10,5) 191 

r 192 

r READ ISP, SURFACE PARAMETER 193 

r lSPa2 UPPER surface calculation 199 

C ISP*0 Lf’^'FR SURFACF CALCULATION 195 

C 196 

RFAOCI.IOO) ISP 197 

MRITEC2, lot ) ISP 198 

iOJ FORMaK;/* TSPaAl?) 199 

IFdSP.EO.a) GO TCr 107 150 

no 105 Mat , INTO 151 

IbIMTO.K)+1 152 

XA(I 5*-xahldCN) 155 

105 ZA (Da-ZAHLOfM) 159 

GO TO 110 155 

107 DO 109 NBl,lMT0 158 

XA(N)aXAHLC'N) 157 

109 ZA (N)sZAhLD(M) 158 

1 10 cnwTIMUE 159 

160 

r CALCULATE DISTANCES AROUND SURFACE CORRE8PONOING TO DISTANCES 181 

C ALONG CHORD LINE, 182 

CALL GEOMTRY CISY,RHO,CH) 183 

WRITE (2, 90 3) tN,X A(N) ,ZA(M) ,8TH(N) ,TH(M) ,FSTM(M) ,FZTH ( N) , MbNI , N2 ) 169 

903 FORMATnHO,9XlHN,5X5HXA(N) ,11X5HZA(N),11X8H8TH(N) , lOXSHTHtN), 165 

110X7HF8TH(N)reX7HFZTH(N)/ 166 

2(lX,I5,6E16,e) ) 167 

168 

C REAn TMO OUTPUT PARAMETERS, NO. OF POINTS IN OUTPUT LIST (IF IFPTPl) 169 

C and output interval in X (IF IFPT»9), 170 

C IFPTpI. full print out every step, IFPT»2', full P/0 at points in list, 171 

r IFPTpS, FUU p/0 AT VFLOCTTY DATA POINTS, IPPTP9, FULLP/0 AT EVERY 172 

C NTH STANDARD 0/P STEP,DX, IN Z DIRECTION EVERY INTI POINT 18 PRINTED, 173 

179 

75 READ(1,106)IFPT,TNT1,MLIST,DX 175 

NLlSTPlBNLISTtl 176 

106 FORMAT(3I5,2F10,5) 177 

wRITE(2,909yiFPT,lNTl,NUST,DX 178 

909 FORMAT(1H0,5HIFPT«,I5,2X,5HINT1«,15,2X,6HNLI8T«,I5,2K,5HDXb,F10,5) 179 

180 
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OBIGINAO PAGE IS 
0?! POOB QUALITY 



IFnFPT,Nt,2)G0 TO 310 181 

C «FAD LIST OF OUTPUT POINTS 182 

read (1,501 HnPX(M),NB2,NLl8TPn 183 

501 FDPMATCeFlO.S) IBU 

WRITt(2,a05HOPXCN),N«2,NU6TPn 16? 

a05 FnRMAT(lHO,6HOPX(N)/flH , 6 (FI 0,5»2Xn ) 186 

310 CONTINUE 187 

186 

C RFAO 8UR"TITLE 16R 

73 PEAOn,lSl)C 150 

WRITE(2,«2nc 191 

421 FORmAT(1HO,3HCd »RA103 192 

195 

f REAO angle of shear (IN DEGREES! ANO TRIPWIRE OIAHETER, 19« 

195 

REAO(l,50nP8I,DTRIP 196 

wRlTEC2,a063P91 ,DTRIP 197 

i )06 FnniHAT{lHO,«HPSIs,F10.5,2X,6HOTRIP«,Flo’.51 198 

199 

PSIBDTPaPSI 200 

201 

C0SP*CO8 (PSD 202 

StNPw8IN(PSI) 203 

204 

C READ NUMBER OF VALUES OF PEYNOLOS NUMBER ANO nEFINlTION PARAHETER 205 

C IFRbi'. RN*0L*SEC (PSD /Nil, 1FR»2, RNnKIL/NU, IFRb3, RNbOLPCOS (PSD /NU 206 

207 

REAO( 1 , 1 OOHnTRL, IFR 208 

REA0(i,501)(RNL(I),I«.1,TNTRU 209 

wRITEC2,407HNTRL,IFR 210 

4107 FORMAT{lHO,6WINTRLP,I2,2X,aHIFRB,12) 211 

wriTE(?.,40R)(RNLCD,I*1#INTRL) 212 

408 FORNATnHO,6HRNLrn/(lK , 8 (F 1 0 . 1 , 2X ) D 213 

214 

c read velocity data Parameter and no. of pairs of Values to follow 215 

c TNTVsl, II FOLLOW, INTV*2. UwSEC(PSn FOLLOW, TNTV«5. CPS FOLLOW, 216 

217 

98 REAOn ,-100! INTV, L 216 

LPlsL+1 219 

READ(l,5iOHXvn),UH(n,IR2»LPn 220 

IF(TSP,E(3.2) GO TO 112 221 

DO 108 Is2,LPl 222 

lO'a XV(I)«^XV(I) 223 

112 ICOU.nTbO 224 

wRITE(2,409)1NTV,L 225 

4109 FORMAT C1H0,5HINTVi»,I2,2X,2HL«»ISJ 226 

IF(INC,EQ.O) PRINT 5q3 227 

503 F0RMAT(/* FLOW IS COMPRESSIBLE AND THE FOLLOWING UM(1) 18 THE MACH 228 

1 NO, DISTRIPUTTON NORMAL TO THE leading EDGE (MENJw) 229 

WRITE(2,410)n,XVtD,UM(l),lR2,LPl) 230 

IF (INC, EO,n GO TO 505 231 

AMINFRbI ,/AM lNFSD 232 

DO 504 I«2,LP1 235 

504 UMdJwUMCnwAMlNFR 234 

505 CONTINUE 235 

410 FORHAT(lHO,4X!HI,3X5HXV(nf5X5HUM(I)/nx,I5,2FlO,5n 236 

237 

UM( l)»0,0 . 238 

239 

C CONVERT velocity DATA TO U. , 240 

30 : 



ortorun r* n or>oor»r»f* 


IFCIN'TV.MC.nCALU VELOCtS (INTV.COSP) 2PJ 

242 

C RtAD VELOCITY GWaDIE^jT AT ATTACHMENT LINE AND, IF VG«AD > 0 AND 243 

C INT3 NOT a 0 , DISTANCE FROM LFADIMR EDGE TO ATTACHMENT LINE • 24<l 

IF (INC. EG. n GO TO 422 245 

246 

FLOR T8 cnMPRjrSSlBLE, READ IN MGR Ap«DMEN/DSCOMP, THEN COMPUTE VIA THE 247 
STEHARTSCIN transformation THE TRANSFORMED INCOMPRESSIBLE VELOCITY 248 

gradient, VGRAO, where 249 

VGRAnap{UPRlHE/nN) /OSINC 250 

aVGRAD/AMlNF5Di*DSCOMP/D8INC*8EC(PSn 25l 

252 

READn,S0n MGRAD,yVCn 253 

IFCTSP.ED.a) GO TO U3 254 

yvtna-xvcn 255 

113 wRlTE(2,42ft) MGRAO 256 

42fi format (/* MGRADbaE16,R) 257 

RATIOBl ,/(! .♦AMFS3D*C08P**2) 258 

VGRADBMGRAD/AMlMF3D*RA7I0**GAri/'C08P 259 

GO TO 427 260 

422 REAOCl,50n VGRAD,XV(n 261 

427 wRlTE(2,4in VGPAD,XV(11 262 

411 FDRMAT(lH0,6HVGRApB,Fl0,5,2X,6HXV(l 1«,E16.S) 285 

IFCVGRAD.GT.0. 1 GO TO 122 264 

265 

IF(ISy,EQ.O.OR,IMT3,EO.O)CO TO 14 266 

267 

read DATA FOR 8/R TO COMPUTE (OU/OX) A. L. 268 

REAOCl ,50UALPHA,8t,S5 269 

WRITE(2,4i2)ALPHA,Sl ,83 270 

412 FnRMAT(lH0,6HALPHA«,F8.5,2X,3H8lB,FB.5,2X,SHSSB,FS.5I 271 

alphabalpha/cosp 272 

GO TO 13 273 

274 

14 WRITE(2,3125 275 

312 FORHAT(55HOVELOCITY GRADIENT MUST BE SPECIFIED WITH THIS GEOMETRY) 278 

GO TO 21 277 

278 

compute velocity gradient at ATTACHMENT LINE 279 

15 CALL VGRaOATCALPHA,RHO,81,8S,VGRAD,XV(1)) 280 

XV(l)Byv(l)*CH 281 

282 

INTERPOLATE FROM TABLES TO FIND DISTANCES FROM ATTACHMENT LINE TO 283 

POINTS at which velocity DATA 18 GIVEN, 284 

122 CALL 8THFRHX (L^H,XV,SXV,THXV,D8DT,INT3.CH,XVCn,SATTi 285 

PRINT 5U,8ATT 286 

5U FORMATt/* non-dimensional DISTANCE FROM LOWER SURFACE TRAILING EDO 287 

le (IF UPPER SURFACE IS TO BE COMPUTED, I8Pp2) , OR FROM UPPER SURFAC 288 

2EA7W TRAILING EDGE CIF LOHIP SURFACE IS TO BE COMPUTED, ISPpo) TO 289 
3ATTACHMENT LINE a 8ATT/CM •*Ei6‘,8/> 290 

IF (INCiEQ.D GO TO 508 291 

292. 

FLOW IS COMPRE8SIPLE, USE THE 8TEHARTS0N TRANSFORMATION IN 293 

8/R 81NCFRS TO TRANSFORM THE ACTUAL SURFACE DISTANCE, S, INTO AN 294 

EOUIVALENT incompressible DISTANCE, 81NC *‘>5 

296 

CALL SINCFR8(LP1,8XV,8XVINC,UH,0SCD8I1,DSCDSI2) 2R7 

GO TO 507 298 

508 DO 506 N«1,LP1 299 

506 SXV1NC(N)»8XV(N) 300 




ri o r> o o 


r,D TO 509 501 

502 

FMOM SXVIMC and SXV(SXVINC) COMPUTE FSVSIWC FOR INTERPOLATION 305 

FROM INCOMPRESSiBLEf TRANSFORMED (8TEWART80N) PLANE TO 300 

compressible* physical plane 305 

306 

507 CALL CSG(SXVINC*SyV*FSVSlNC*LPlfD8C08n, DSCDSI2) 307 

509 CONTINUE 308 

309 

C '“BITE main and sub titles 310 

WR1TEC2,83)B,C 3ll 

63 FO«HAT(1H1,8AIO/1HO*6AIO) 312 

313 

wHlTE(2,oitt) 3ia 

liiu FORHATflHO.SlHREYMOLOS NUMBER DEFINED BY RN» ) 315 

IFnFR,E«.nHRlTEt2f015) 316 

IF(IFR,EQ.2)NPITE(2*R16) 317 

IF(IFR,EO,5)HRITE{2*017) 318 

015 E0RMAT(lH+,51X,l«HQL*8EC(PSn/NU) 3l9 

016 FORMAT(lH + i3lX,lCiHBL/NU ) 320 

417 FDRmatC1H4,31X,10HOL*COv (P81)/NU) 321 

91 wRlTE(g,92) VGRAO 322 

9? FORMaT(39HOVEL.OCITY GRAHIENT AT ATTACHMENT LINE* ,F8,2) 323 

324 

L*L+1 325 

326 

LPl*L+l 327 

C FROM theta and UCTHETA) COMPUTE FUTH FOR VEl.tlClTY INTERPOLATIONS, 328 

CALL CSG(THXV,UM,FUTH.L.n80TAVGPA0AC0«P, (i 'M (L)-UM f |,-l U/ (THX V (U- 329 

ITHXVCL-nn 330 

C 331 

r WRITE TABLE OF VELOCITY DATA 532 

CALL X8CPPNTCINTV) 533 

334 

nl1ST*nlIST+1 335 

336 

C COMPILE LIST OF OUTPUT POINTS IF REQUIRED, 337 

IFdFPT.NE.l )CALL PLI ST { IFPT, NL! 8T * IXV , 8XVINC ,L » D8 , TNT3 * CH , XV ( 13 , 538 

IPX, INC) 339 

340 

C IvS NING SHEARED 541 

IF(PSI,Lt..OOOOnGO TO 26 342 

345 

344 

C GUESS AT velocities AT ATTACHMENT LINE 345 

26 AbJ 346 

DO 9 N«2,J 347 

UMl (N)*(N.I )/a 348 

9 CONTINUE 349 

350 

C CONSTANT 351 

AO*»0,5/(DZ*OZ3 352 

553 

C SET STEP COUNT 354 

M8l , ,355 

MS*1 . 356 

357 

IFCONBO 358 

359 

r boumdary cdmoitions including suction up injection, 360 


32 



UMl n)«UMC(n«UH2M )«VM2(1 >bO, 361 

WCOKilrn , + AMf 83nikf:oSP<t*25'**(-.S)/(l .4AMFSSP1 562 

'*'M(n»HWAtL(n/SCRTcvuf»AD*kCDSP)*WCCIMl 363 

UW2tJ + n*VH2(J+nBj , 364 

365 

C GO TC B/P wUV TO COMPUTE W, tl AND V AT ATTACHMENT LINE, 366 

367 

CAtL WUV(1,SC2)) 368 

369 

AN0LF1«ANGLE2*,5*PI . 370 

SCQMPsO. 371 

AME3n«$RRTC(AMlNE3Diir*2*SlNP**i2)/(l , + AMF83D<»C08P**2)) 372 

373 

C WRITE ATTACHMENT LINE PROFILES 574 

375 

C leading edge contamination test 376 

377 

»*-RITE(?,t0<O 378 

104 F0RMAT(lHl/iy65(2H*«)//6aX6HClLl7PUT//lX65C2H*-n 579 

caul CONTAM(VRBAD»07RIPiCH,C08P,STNP,THETAl,RAT!0) 380 

call PRlNTtXV{n»O,,0.,VGRAn*CO8P,J»D7,T.NTl,l ,0,,0,ANGLE2) 381 

362 

DO 16 N*1 , INTRL 363 

KMAX(n)bO,0 384 

18 CONTINUE 385 

8(1)B0, 386 

SCOMPlaiO, 387 

Um«UMcn 386 

DSZbDS 3BR 

390 

NEXTBl 391 

SNEXTb«1,0 392 

INTHOLDbO 393 

LAST«0 394 

ISTbI 395 

CFXIbO, 396 

COFXbO, 397 

CDFXINFbO, 396 

CFVlB2,-*SINPAr(VGRAD*cnSP/RNL(l))*A,5*DVDZ 399 

UOO 

C ADVANCE STEP COUNT 401 

1 1 MbM + 1 402 

LCbO 403 

404 

r CALCULATE LENGTH OF NEXT STEP IN X DIRECTION 405 

406 

10 CALL STPLNTHCSNEXT, 8fINTH0LD»D8,08l,D82, NEXT, NLI8T>LAST» 407 

1IFPT,LC,ILP»U,USTEP.DU08,X>8H»ITC,ANGLE2»WF) 408 

409 

c IS COMPUTATION TO END ttlO 

IF CLAST. £0,2)60 TO 239 4H 

412 

C DDES THIS STEP END AT A LISTED OUTPUT POINT «13 

IFCILP.NE.DGO to 49 414 

415 

C STORE UCM-2), U(M-n AND V(M« 1 ) PROFILES WHILE A STEP ENDING AT A 416 

C LISTED OUTPUT POINT IS COMPUTED 417 

on 46 Nb2,J 418 

UM1H(N)bUM1 (N) 419 

UM2H(M)bUM 2(N) a»0 



U(> CONTINUE 

continue 

C VARIABLES IMOEPENOENT OF N 

UAVG»WF*U(2)4-M ,-WF)«!U(l ) 

G«WF**2/4.* (SHy (UAVG*OSl )*(U(2)-U(n )-l.) 

CONl sl ,/(:iZ*r*2 
CON2*n ,-WF)/WF 
CDN3«8H/(DS1*UAVG) 

CONiaaM ,^WF)*8H/DS1 

COMSsf n ,-!^F)ywF)**2 
CnNfeBSH*nZ/(2 .Ti»>fF*DSJ *UAVG) 

C0N7»SHy081 

c variables dependent on n 
DO 5 N«2,J 

Gl«G*(UM2CN + n-UN2(N-} n 

AND(N)«C0N2*Gii.(N-1)fUM2(m) 

WNOCN)»CONfe*C )*U(2)A2.*n,-WF)*U(l ) )*(UH?(N)AUM|(N-n ) 

.•WFlikDZy f2,*.^F**3T*r,* CUM2(NT-UM2(N ■^1 n 

IFt'.lF.UT, ,80n^NCi(N)a.S*WNO(N3 

RN0(N5*-4,*WF<!A0 + C0N2*(N-n*G*CUH2(NAn»UH2(N«l ) ) 

IF(WF,GE,.75) GO TO 512 

BND(N3 ■CON3*WF* t3,-a,*WF)*U(2) AUM2(N5+BNDfN) 

512 DN0(Nl*C0N3*CU(2)-uai+((2»*<^F»l .T^d'.-k^FUUtllAa'.FCl.-wnAAaPLKn 
i)*uM2(N3*uM2CNn4C0NiAa ,-iNF)*(UFa(Ntn-2.*uM2(N)+uN2(N-n ) 

2-C0N5* (N-i ) AUH2 tN}F(UM2CN+n-UH2 (N-n 3*G 

BDOtN)aCDN4#UM2(N)'F5» *wF*CON1 

ODO(N)»CON4AVM2(S5*UM2(N)+Cl,-wF)*CONr*(VM2(N+n-?,*VM2(Nl+VM2(N-l 

n )"CON5*GA(N-n*UH2(N)*(VH2(N+13-VM2(N-ni 

C EXTPAPOLATE TO ESTIHATE UM PROFILE 
IF{N,GT.3)G0 to 2 
UPG(N)*UM2(N) 

GO TO 3 

2 UHG(N)aUH2(Nl«0Sl /OS2*(UH2(NJ-UMl (N) ) 

3 CONTINUE 

C complete THE SPECIFICATION OF SURFACE SUCTION OR INJECTION VELOCITY, w, 
IFC18LC.EO.03GO TO 201 

NC0N2aSQRTa,F.5*(GAMMA-l,)*AME30*F2)y (1 ,AAMF830)**1 .5 
8HC0MP«WF*SC0MP+tl ,-WF)*SCOMPl 
TFtSPSCMS) .GT.SHCOHP) GO TO 200 
MSaHS+l 

200 WM(t)»iwWALLtMS3*(SH/UAVG)**.5*(NFAWCON8+(l,*WF)*WCOHn 
GO to 202 

201 WM(i)pO, 

202 CONTINUE 

C COMPUTE WM, UM AND VM PROFILES 
CALL WUV(M,8(2)) 

f. COMPUTE SKIN FRICTION DRAS, BASED ON REYNOI PS NUMBER DEFINITION, IFR»2 
C AND THE FIRST VALUE IN REYNOLDS NUMBER ARRAY, RNLCIl, 


021 

422 

423 

424 

425 

426 

427 
426 

429 

430 

431 

432 
435 

434 

435 
4 36 
437 
458 

439 

440 

441 

442 

443 

44 0 

045 

o4ti 

447 

04S 

44Q 

050 

451 

452 

453 

454 

455 

456 

457 

458 

459 

460 

461 
4 62 

463 

464 

465 

466 

467 
46S 

469 

470 

471 

472 
475 

474 

475 

476 

477 
47s 
479 
460 


3h 



c CFX«TAUK/{,5*»H0INF*raiNF*<t2) «5l 

C CFY*TAUy/C,5*RH0INFiM3INF**E) 

C CDFX ■ 8KTN FRICTION DRAG COEFFICIENT IN DIRECTION PERPENDICULAR TO LE ^03 
C BASED ON CHORD HEA8URED PERPENDICULAR TO LE 

C CDFXINF • TOTAL SKIN FRICTION DRAG COEFFICIENT IN FREESTRFAM DIRECTION <105 
C BASED ON CHORD MEASURED PARALLEL TO fREEBTREAM <J6fe 

R67 

TFCILP.EO.DGn TP 203 

IC0UNT«1CDUNT+1 R09 

D8COHP«8COMP-SCOHP1 . RRO 

CFX2«2,ArU(2)*CU(2I/tRNL(nAS(2)n**t5*DUD'Z <<91 

CFY2»2.*8INP<,tlJ(2)/CRNLCn*8{2n)A*,5*nVDZ «9Z 

DCDFX«,5*08COHP*(C08(ANGLEn*CFyi+C08(ANGLE2)ACFX2) R93 

VCFy2CICOUNT)«l00O,*CFY2 U9« 

VCPX2CICOUNn«l OOO.ACFX2 R95 

H8COMP(ICOUNT)«3COKP U96 

CDFX«COFXtDCDFX «97 

CnFXINE«CDFXINEA0C0Ey*CO8P*,5*8TNP*O8COMP<iCCFY2+CEYU R90 

203 CONTINUE «99 

500 

r HAS ITERATION CONVERGED 501 

JFdFCON.EO.ncO TD 239 502 

503 

C HAS STEP length BEEN HALVED TRICE ALREADY 50a 

IF(LC,e0.2)G0 TO 23B 505 

50fc 

C HALVE STEP LENGTH AND TRY AGAIN 507 

IC*LC+1 508 

OD TO 1 0 509 

510 

23B IFCONaJ 511 

LA8T«2 512 

513 

C DETERMINE IF PRINT OUT 18 TO 5E COMPLETE# PARTIAL OR SKIPPED 5l« 

239 CALL IFPR1NT{IFPT#ILP.LA8T»JACKP0T) 515 

510 

IF( JACKPOT. EO.OIGO TO 5 517 

510 

C COMPUTE CROSS-FLOW VELOCITIES 519 

lF{PSl,GT.,OOOn CALL CR88FLW ( J ,OE# U (2 ) ) 520 

521 

C WRITE BOUMDARV-LAVER CALCULATION RESULTS A8 REQUIRED 522 

CALL PRINTex,S(25,UC23#0UD8,J#02»lNTl»JACKP0T,PSl,LC,ANGLE2) 523 

52« 

C CALCULATE. OIMENSIONALISING factors 525 

CALL DMN8I0NeS(2)iU(2)#C08P,DELtAl,THETAl,ILP.CFX2,CFy2»CDFX, 526 

ICOFXINF) 527 

528 

C RE-LAMINARISATION TEST 529 

IFCKMAX(INTRL) .GT’,-0,5)CALL RELAM(U»8i8INP»C08P|KMAX) 530 

C 18 PRINT-OUT complete 532 

IFCJACKP0T,E8,2)G0 TO 5 533 

534 

C CROSS-FLOW instability TEST 555 

IF(P81,GT.,000n 536 

1 CALL lN8TAB(C08P»8t2).U(2),CH»RAT10,AMES0»AMlNF3D5 53T 

530 

C VISCOUS instability test 539 

50 CALL TRAN8(8(2),DU08,THETAl,Ut2)#U(l),IST»JACKPQTV 5«0 
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C IS STEP JUST ENDED THf LAST ONE REQUIRED 
IE(L*ST,EQ,n GO TO 27 

C DDES LAST STEP END AT A LISTED OUTPUT POINT 
lEdLP.EO.nGO TO «7 

umsuia) 

smasiaj 

SCOHP1«08COMP 

HCONiawiCONH 
DSZ*DSl 
CFXIBCEX2 
CEYlaCPV2 
angle UANOLE2 
GO TO 11 

C REPLACE U{M-n AND V(M.l) PROFILES -'<ITN THOSE STORED AT START 

c OF Last step 

U7 DO <ja NS2,J 
u^l (N)PUMIH(N) 

UM2(N)KUMaHtN) 

VN2fN)*VH2H(N) 
ttS CONTINUE 
GO TO 11 

27 i»(RIT£t2#lR) 

IR FORMAKSSHOLAMIMAR flow CALCULATED TO END OF DATA OR LAST POINT HE 
IQUESTED) 

GO TO 21 
5 wRITE(2,20) 

20 FORMAT ( UHOSEPARATION) 

0 READ CUE TO READ MORE DATA OR TO FINISH 

21 READ( 1, lOOIKUE 
hr1TE(2,«13)KUE 

U 13 F 0 HMAT( 1 Hl,«HKUEa,I 2 ) 

C 

r PLOT IM3TRUCTI0NS FOR X ANo Y SKIN FRICTION DISTRIBUTIONS 

C 

MPTSPiCnUNT 

HGTb, 1“F HOTia.US HGT2«,068 HGTS«,OS 

NPloNPTS+1 

NPa»NPT8T2 

YORG«0,S XORGaO, 

XSCALEP.a 

XPSbS. 

YP0*6. 

XDVpO,® XTIC«-1. 

YDVpO.S YTlC*-»i, 

ORGsOt 

CALL 83CALECVCFX2,YPG,NPTT(1,Ii1,p« 1»0RG) 
VSCALEaVCFy2(Wp2)»vCFX2(NP2) 

VCFX2(NPTS + nRVCFYa(NPTS + l ItYORO 
VCFX2CNPT8+a)KVCFY2(NPT8>?)RYSCALe 
HSCOMP CNPT8+13»XOHG 
H8C0MPCNPTS+21BXSCALE 
CALL CALPLTC2 '.,S.p«55 

CALL AXES(0,040.0i90,,YPOrVCFX2(NPl)rVCFX2(NP?)»YTlG,YDV,lH ,HGTl^ 


5AII 

542 

5«5 

544 

545 
mh 
547 
54fi 

549 

550 

551 

552 
555 

554 

555 
55S 

557 

558 
55R 

560 

561 
§62 
561 

564 

565 

566 

567 

568 

569 

570 

571 

572 
571 

574 

575 

576 

577 

578 

579 

580 

581 

582 
581 

584 

585 

586 

587 

588 

589 
890 

591 

592 
595 
59 a 
595 
196 

597 

598 

599 

600 
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i+n fcoi 

CAUL NnTATE(-.5»2,5fH6T»<?H CFXa 1 000 , 90 , , «» 5 602 

CALL NriTATEC-.75,S.5,HGT,9H CPY*tOOO,90'. ,9) 603 

CALL AyR8(0,0,0.0, 0,0»XPG,HSCOmP(NPJ ) , HSCCIMP (NP2 ) , XT IC # XOV , 609 

116H8URFACE D 1ST AfJCE » HGT I , - 1 6 ) 605 

CALL LINE(HSCnMP,VCFX2,NPTS,l,+l,3,,06) 606 

CALL LIMECHSC0MP,VrPY2,NPT8, 1 ,tl .05) 60? 

CALL KiOTATFC.5. YPG,MGT,27HSXIN FRICTinN DT 8-TR IBUT IONS , 0 , , 27 ) 608 

HGLi.ia 609 

XPLTB.5+.5AHGL 610 

YPLTPYPG»,aO+.5*HGL 611 

CALL Mn7ATECXPLT,YPLT,HGL,3.0.,.n 612 

CALL NOTATECl .0,YPG-,00»HOT,3HCFX,0.,3) 613 

YPLTbYPG-,70 + .5*(HGL-.O2) 6U 

CALL WOTATE(XPLT,YPLT,HGL-.02»«,0,,-n 615 

CALL NOTATEC1.0,yPO».70,MGT,3HCFY,0,,3) 616 

P8I0*P81A57, 29577951 617 

EtJCODEC80,a23,LAflFLCn) AMI NF3D , P8IO, RNL ( 1) 618 

425 FORMAT(0HAMlNF3D«F5,3.3X,UHP8I«P4,l,3X7HPMLn)BlPE9,2,i7X) 619 

CALL MOTATE(0.,-l,oo,HGT,LAflEL(niO*rPO) 620 

EMCOnE(80,a26,LABELn n COF X , CDFXI NP 621 

426 format ( 5HCDFXK,F7.5,5X,8HCnFXIWFaF7.5,4RX) 622 

CALI MOTATECO,,-J,25,HGT,LABELCn»0,,«0) 623 

CALL MnTATF(O.,-l.50,HOT,LAfl(l).0.,fiO) 624 

IF(KUEiFQ.nGO TO ?.4 625 

iF(KUE.E0.2)Gn TO 74 626 

IF(KUF..F0.33g0 to 76 627 

IF(KUE,EQ,4)G0 to 75 628 

IF(KUE.E0.5)G0 TO 73 629 

1F(KUE.FQ.63GD TO 96 6J0 

GALL CALPLT(0.,0.,P99) 631 

STOP 0101 632 

END • 633 


SUBROUTINE SINCPRS(J,SXV,SXV1NC,UM,D8CDSI1»D8C0S12) 634 

DIMENSION SXV(l),SXVINC(n»UM(n 635 

COMMON/CO^PRES/INC, AHINFSDi AME30,8INP,Cn8P,GAMMA,QAMl ,GAM2»GAM3, 636 

1GAM4.AMFS3D, RATIO, SCOMP,ZCOMP(170),7(170)»RHOO(170) 637 

C 638 

C USE STEWARTSON transformation TO CONVERT 8XV,THE PHYSICAL OISTAN'CE 639 

0 AROUND THE AIRFOIL AT WHICH MACH NO, 18 GIVEN, TO 8XVINC, AN 640 

C EOUIVALENT incompressible COORDINATE 641 

C 642 

GMliGAMMA«t , 693 

DENOMRbI , /C1,+AMFS30*C08P**2) 644 

SXVlNC(n*0. 645 

CONbCI ,aAMF83D)a*RaM1 646 

FACl»l,tAHF830 647 

DO 10 N«2,J 648 

AME3DN»AMINF3DA*2A(UH(N)»*2»FACUSIMP**2)*0EN0MR 649 

AHE50NlBAMINF3D**2*(UM(N«n Aw2*rACUaiNP**2)*DEN0MR 650 

AMENbI ,+,5*GM1*AME3DN 651 

AMENl»l ,t.5*GMlAAME3PNl 652 

10 8XVINC(N)B,5*C0N*(AMENA*(-GAHiyAAMENl*A(-6AMn )* 653 

1 CSXV(N)-SXV(N»l))+8XVlNCCN-n 654 

AHE3ni«AMlNF30AA2*(UM(n6*2*FACl+8INP**23*nEN0MR 655 

AHE30JBAMINF3nA*2*(UHtJ)**2*FAClf8INP**2) ADENOMR 656 

AMEl«l,+,5*OHl*AME3Dl 657 
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AMEJBl , + .S*GMl*A^'ie5DJ 650 
08CDSIl«AMEi »*Gam 1/C(1M 659 
nSrDSIg«AMEJ**rGAMl/cnK( 660 
RETURN 661 
ENO 66? 


665 

subroutine WijV(M,S) 66U 

665 

C COMPUTES w, u and V PROFILES# 666 

667 

CnMMON/SURWUV/J,ITS,TnL.lFCON,OS,nSl ,nZ,wF,8H,l.lHGH,A0,C0N3,C0N6» 668 

1C0N7,G, ANt; 170 5 .ROOtnO) , RN ( 1 7 0 5 * BNO (1 7 0 1 , RND ( I 70 ) , CN (1 70 ) , 000 (1 70 ) 669 

2,DNn 705 ,BCn,DNOn7O),nNnn70),Ul^l(UO).UM6n7Q),WM(1705»WKl(17n), 670 

5rthn(l7n) ,U (0) , ANn(l70) . 671 

CnPMnN/RESlJLTS/WM2 ( 1 70 5 ,UM2( 1 70 5 ,VB3( 17m ,0ELTA1* THETA 1»N0,PUDZ 67? 

l.nvnz 673 

COMPOfJ/cnMPRESyiNC, AMINF5D»AME50»SINP,c08P,r,AMHA,fiAMi,GAM2,GAM3f 670 

lGANa,A^^Fv'55•’3, RATIO, SCOMP,ZCC!Mp{l70),T(17(m,«H0na7n) 675 

DIMENSION UP(i7Q5#yC17Q> 676 

677 

Nn«o 670 

67R 

C evaluation of u. 680 

601 

11 DO 2 N«a,J 602 

IF(N,NE.nGO TO 1 60S 

r, w at attachment line 68a 

WM(N5«WH(N-l)»0.5*OZi*(UMi{N)tUHl(N«in 605 

GO TO 2 606 

667 

C W AT X(M-i/23 60S 

1 FACTs-WF*cnN6<ii(2,*U(25 + Cl , -2 ,*WF) /wF*U ( 15 5 * (UHg ( N) fUMG ( N. 1 3 ) 669 

l-(a*N-65*.5*DZ/WF**2*(UHG(N3-UMR(N.n5*G 690 

IK(WF.LT,.501 3 GO TO 50 691 

WM(N3BWM(N-l)-(l,-WF3/wF*(WMl(N3-WMl(N«n5+HNn{N)A.FACT 693 

GO TO 2 695 

50 WM(N3«WM(N-1 3+HNO(N57.5FFACT 69« 

2 CONTINUE 695 

696 

C evaluation of LIMiN 697 

698 

DO N«a,J 699 

IFCM.NE.DGO to 3 700 

701 

AN(M3«0,S«WM(N3yDZ72*AD 702 

0N(N3»UMl tN)-a*AO 703 

CN(N3B-AN<N5+a*AQ 70a 

DN(N3«1 705 

GO TO a 706 

707 

3 IF{WF,LT,.50n GO TO 51 7Q8 

AN(N3«2. *WF*A04'.5*WF/nZ*(WF«(WM(N3tn ,-WE3*WMl fNn + fN«l)*6tUMGCN) 709 

1+AN0CN3 710 

nN(M3«DNOfN)+CON3*u( l )Nk(l .-WF3*C«’,*WF-1’. J*UH2 (N 3 *UM0 (N3*«5/DZF ( 1 , - 711 

lWF3*CUM2(Nfl3-UM3(N-i53*(WF*WM(N)+(l,-WF)*WMl(N)5 712 

GO TO 52 713 

5) AN(N3« An4..25/DZ*WM(N3 + tN-n#G*UMO(N5+ANO(M) 7H 
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l5M(Ki)*Dt»0(N!U.5*Cntvi3*IJ( n*UM2(ig)*U(*G(M)-,85/DZ*(UM2(N + n»UM2(N-l 3 3 7 15 

t*i'<M(N3 716 

5? RN(M)*CP''i5*t2,*WF**?*U(2)+wF*(l‘.-2,*y».p)*Ucn)*UHG(N)+BNOCK‘3 717 

CN(M)«-ANi(N) + a,*KiF*FO 718 

FACT«»CON3*WF*(5,-«,*MF)*IJC23*UH2(Im3i>UM[;(W3 719 

TF(WF.LT,.75) Gn TO H 720 

D^l(N)liD^(^13+PiCT 721 

U CONYINUF 722 

723 

LP(J3bBN(J3 72« 

00 5 K«3fJ 725 

NhJpK+5 726 

UP(Ki«l)BflN(N«n-AW(M»\ )*CNCM)/UP(M3 727 

5 CnvTJi^UF 728 

729 

Y{J)BDNtJ3-AM(J) 730 

DO 6 Ks!3,J 731 

NB,r«K+3 732 

Y£N|-l3*0M(N-l)-A>J(»J"n*Y(Si)/UP(N) 735 

h CONTINUE 730 

735 

UMO(a3»Y(23/UP(2) 756 

DO 7 Nfi3»J 737 

UNSCN3»(V(m3«cn(N3*U^'RCN-1)3/UP(N) 738 

7 CONTINUE 759 

7*l0 

C COUNT number of iterations 7«1 

nObMO+1 742 

TFtNO.LT.3) GO TO 22 745 

7ua 

lF(T0L,GT.AeS(UM6H-UMr.(:2n/0Z)GrT TO 8 7«5 

f CHECK number of iterations 7A6 

iPtNO.GE. IT83GO TO 12 7«7 

7^8 

C STORE U NEAREST 8URFACF FOR CONVERGENCE CHECK 7<i9 

22 LIMGH*UMG(2) 750 

JFiM.NE.nGO TO n 751 

752 

DO 37 N»2, J 753 

UMt (N)eiUMG CN) 75« 

37 CONTINUE 755 

GO TO 11 756 

757 

C ITERATION HAS NOT CONVERGED 758 

12 IFC0NB2 759 

RETURN 780 

761 

C iteration has converged 762 

ft IFCONbI 763 

76^J 

C EVALUATION OF VH,N 765 

786 

00 16 Nb2, J 767 

768 

lF(M,GT.nGO TO 21 769 

0NP(N)B-a*AO 770 

DND(N)*0 771 

GO TO 16 772 

773 

21 PND(N)«POOCN3+WFACriN7*UMGi:N3 
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IFtwp.LT, .son GO TO 53 . 775 

FACT«,5/DZi»(l.-wn*(WF*WH(M) + U’.-v<JF)*!WMi (Nn 77fc 

GO TO 5« 777 

53 FACTb,25/DZ*wM(N) 77» 

5/4 DNO(M)«nnO(W)+CnM7*WF*Vt^2tN)*UHr,(M)-(VH2(KJ+n-VM2(N"n )*(PACTTC1 ,• TV* 

lK,F)/WF*(W-n*G*UMR(N) ) 780 

78i 

COMTINUE 782 

7 83 

UP(J)*BKiO(J) 76a 

DO 17 Ka3,J 765 

NaJ-K+3 786 

upt‘4-nioHNncK.-n-AM(N(-n*c'''(N)/up(fJi 78? 

17 continue 788 

789 

Y(J)«0N0(J3-4N(J) 790 

no 18 Kb3,J 791 

NaJ-K+S 792 

Y(N-naDND(M-l)-AN(N«n*Y(N)/UP(N) 793 

16 continue 79a 

795 

VM2(2JaY(23/UP(25 796 

00 19 Na3fJ 797 

VK2tMBCYtNl«CN(N)*VN2CN-in/UPfNJ 798 

19 CONTINUE 799 

80 0 

00 23 Nan J 801 

C STORE W PPOFILE FOR POSSIBLE PRINTIN 5 OUT, 802 

^1^/^(:H3aWM(^45 805 

0 LI PROFILE AT ENO OF STEP BECOMES PROFILE AT START Of NEXT STEP, 806 

UMl (N) hUM2(N) 805 

UM2(N)bUmgCN) 806 

wMl (NlaWM(W) 807 

23 CONTINUE 808 

C 809 

C COMPUTE TEMPERATURE RATIO. T(N)aT/TlNF AND DENSITY RATIO, BIO 

HHODCNlaRHO/RHOlNF, AND THE COMPRESSIBLE (ACTUALJ NORMAL 811 

COORDINATE, ZCOMPtN)aZ/C 818 

813 

IFtM.GT.nGO TO 25 «ia 

U(25a0, 815 

RATlOal ./(1 ,yAMFS3D*C08P*a21 816 

25 ZCONaRATIO**GAMl 8l7 

DELC0NaRATI0<i*GAM2 818 

RH0CONpRAT1O**(-GAM3) 819 

CFCONXaRATlO*A(-GAMa) 820 

CFCONYaRATIO**(-GAMl) 82l 

UCONbRATIO>*(-,53*U{21 822 

ZCOHPtnaO, 823 

jPlBj+1 82a 

00 26 Nal.JPl 625 

URATIObUCON*UH2(MJ «28 

VRAT10»SINP*VM2(N) 827 

T(N3«I1.+AMFS30*(1,-URATI0aa2-VRAT1O**2J «28 

RHOO(N)pHHOCON/T(N) 829' 

IF (N, £0,1) GO TO 26 830 

2CaMPCN)«ZC0N*.5*(TCN) + T(N-in*DZ + ZC0MP(N-n 831 

26 CONTINUE 832 

833 

C CALCULATE DISPLACEMENT ANp MOMENTUM THICKNESSES, 



DF.LTAlBl.O 855 

THRTAImO 83fc 

WPI 837 

IS«1 838 

«0 M«N*1 839 

DELTAl»DELTAl+(l-UMg(M) )*(3+IS) 840 

THETAl»THETAl + (UM?CfO*(l-UH?(N)UC5tI8)) 8«l 

I8S-I8 8<*2 

IFCUH2CK') .Ge.l ,0-TnU*OZ)Gn To /^1 8U3 

IPCN.IT. JAJ )G0 TO «0 . 8<ia 

f5eLTAl*nE!.TAl*DZ/3.0 8ttS 

THETA1pTHETA 1*DZ/3,0 8tife 

ZERGepPLOATCN-n*nZ 8«7 

DEI,TA1,iiZC0^'PCM)-nELCnN<r{ZFDGE-0ELTAn 8«8 

THeTAlnnELCON*THFTAl RUO 

850 

C estimate (DU/DZ) ANO (DV/DZ) at ZpO, 851 

0UDZ*(2. AUM2C2)-0,5<rLlMg(3) )/0Z8CFC0NX 852 

DVnZE(2,0#VMg(2),n,5*Vh2C3))/nZ*CFCnNY 853 

85^ 

RETI.Jf?lM 855 

85fe 

END 857 


858 

FUNCTION THXCX) 859 

880 

C THaNSFORws X TO THETA 881 

882 

ARG«e2.AAB8(X)-l , 863 

IF CARS, LE, 1 ,03 GO TO 2 860 

aRITE(2,3) X,ARG 865 

3 FORMAT (///A ERROR Nfl, 2 OETECTEO BY ACOSIM IN FUNCTION THX*/ 866 

I* Xp*B 16,B,* ARGi8A,Et6.8//) 867 

1F(ARG.GT,1,,AND.ARG,LT,1.0000001) ARGPl, 868 

? THX»AC0-8(ARG) 869 

IFCX.LT.O, )Gn TO 1 870 

THX«6.2831853-THX 871 

1 RETURN 872 

END 873 


B7« 

FUNCTION XTH(THETA) 875 

878 

c transforms theta to X *77 

878 

XTHpO.SaCIaCOSCTHETA)) 8?9 

IF{THETA.GT.3. 1015926)60 TO 1 880 

XTHp-XTH 681 

1 RETURN 882 

ENO 883 



SUBROUTINE C8S(X,Y,F|NGP,|l,Fl) «86 

887 

r. GENEMATFS SECOND DERIVATIVES FOR USF IN C8I (CUBIC SFLlNEJ 888 

889 

dimension X(n.Y(n,F(15,EC388),G(385T 690 

891 

N2aN0P»l 892 

DO « N»2,N2 893 

E(N)a2*(X(N + n-X(H-in 89« 

F(N)s<,*fYeN+n-Y(Nn/(x(N+n-v(Nn-8*cy(N5-y(N-t))/(x(N)-x(N-i)) 895 

il CONTINUE 898 

897 

Nl«? 898 

NSaNOP 899 

900 

E(n«2*(X(2)-X(nj 901 

P(n*8*(Y(25-V(in/(X(2)-X(in-6*F.l 902 

E (N0P)b 2*(X(NDP)»X (NOP^n ) 903 

P (NOP)s«ib* CY(HPP)*Y(NnP»n )/(X(NOF)»XtNOP»nU8*FT 900 

905 

DO 8 NbNI , N 2 90fc 

G(Ml*(X(N)-X(N-m/E(N«n 907 

E(N)»E(M)-GCM)*G(N)*E(N-n 908 

P(NlBP(N)-G(N)*FCN-n 909 

8 continue 910 

911 

F(N2)aF(N2l/E(N2J. 912 

on 9 N3«N1,N2 913 

MbMI+MP-NJ 91 « 

F(N-n«F(M-n/ECN-n-G(N)AF(N) 915 

9 CONTINUE 916 

return 917 

END 918 


919 

SuaPOUT-INE C8I (X,Y,F,NOP,XI,YI ,YX) 920 

921 

C CUBIC SPLINE INTERPOLATION 922 

923 

DIMENSION XCn,Y(n,FCn 920 

925 

DO 12 N*2,NCP 926 

lF(XCN)*Xn 12a2»n 927 

12 CONTINUE 928 

929 

NaNOP 930 

931 

13 Al«0,5AF(N-n*CX(Ni-XI)*(XCN)«XI)/(X(N)-XtN-n) 932 

BlB0,5*F(N)*(Xl-X(N«n)*(XI-X(N-lIl/(X(N).X(N.in 933 

ClaY(M-l)/CX{N5«XfN-m-f (N»n*tX(N5»X(N-l))/6 93« 

Dl«Y(N)/(X(N5-X(N-rn-F(N3*(X(N)-X(N»tn/6 955 

Yla(Al*(X(N)-XI)+Bl*{XI-X(N-n))/3 + ClF(X(N)BXI)>0U(X!-X{N-t)) 936 

YXn-Al+Bl-Cl+Dl 937 

RETURN : 938 

end . 939 


42 



SHIRROUTIMF CONT4M(VGRA(5.nTRIP»CH,CQ8P,SIMP,THRTAl ,RATTO) 


r TEST FOP COWTAMINATIDN AT ATTACHMEMT LINE 
Cni^MON/TEST/RNUnO) ,INTPL, IFR 

Ifl eoRmaTCI X»3RH**a LEADING-EOGE C0''IT AKINATICIN test *A*) 

DO 1 N»l,IMTRl 
WSITE(2#10)RNL(»«) 

c SCALE REYNOLDS NUMBER TO STANDARD FORM, 

C RTHETAaOlNf aSINPaTHETA/NUEDGE 

IECIFR.E0.1)RN»PNL(N5 

IEnFR,E0,2)RN»RNL(N)/C08P 
IF(TER,Et3,3)RN«PNL(N)/C08P**2 
RTHETA«SINP*SQRT(RN/VGRAO)*THeTA!/PATtOA»J .5 

ARITECa.inPTHETA 

IF(RTHFTA-10n.J?,2,J 

? WR1TEC2,12) 

GO TO 1 

3 IF{RTHETA-2U0,)i(.«f5 

4 IF(riTRIPi..OOOOn 

6 DCPIT»CH*47*30RT(RTH€TA)/fRNA8lNP*C08P) 

WRITE(2. 13)DCR1T 
IF(nCRIT-OTPlP)7,l()«r6 

7 f'RITEf?. 15) 

CO TO I 

U WRITE(2.17) 

GO TO 1 

5 WRITE(2.J4) 

1 CONTINUE 

RETURN • 

10 FORMAinSH REYNOLDS NUMBER* ,ElO,0) 

11 E0RMAT(1H+,52X,7HRTHETA",F8.1) . , 

12 E0PMAT(1H+,67X,34HND TURBULENT CONTAMINATION AT A,L.) 

15 FORMAT(lH+,5lX,fcHDCRIT*»F8,<t) 

U F0RMATaH+,67X,31HTURBULENT CONTAMINATION AT A,L.) 

15 FnRMAT(lH+,67Xi36HTURBULENT CONTAMINATION AT TRIP WIRE) 

17 FORMATC 1H+, S7X,40HTURBULENT CONTAMINATION POSSIBLE AT A,L.) 

END 


940 

941 

942 

943 

944 

945 
94(h 

947 

948 
949I 
950! 
9511 
952i 
953‘ 
954' 

955 

956 

957 

958 

959 

960 

961 

962 

963 

964 

965 

966 

967 
966 

969 

970 

971 

972 

973 

974 

975 

976 

977 

978 

979 

980 

981 

982 

983 

984 

985 

986 


SUBROUTINE CR88ELH (J»PZfU) 

C CALCULATES CROSS-ELOW AND STREAH-FLOVi PROFILES AND THICKNESSES, 

C0HM0N/RESULTS/WM2(170) f UM2n70) »VM2n70)iDELTAl ,THETAl ,NO,DUOZ 

I ,DVDZ. 

COMM0N/CPO88V/8VC170),CV(170),80TfCOT,CVM 

COMMON/COMPRE8/INC,AMTNESO» AME3D*8INPiC08P»eAMMA,6AMl,6AM2,GAM3i 
l(;AM4,AMFS3D.RATIO,8COMP,ZC0MP(l70),T(170),RHOOn70) 


987 

988 

989 

990 

991 

992 

993 

994 

995 

996 




997 

C VELOCITY ftT EOGP OF f\DU^iO&«Y LAYER 998 

SVJaSQRT(U**2/RATI0+3TMP**H) 999 

1000 

r. calculate SIMCTHETA) where Theta is angle between flow at edge of 1001 

r 80UN0ARY LAYER AND THE PERPENDICULAR TO THE LEAOINg EDGE. 1002 

SINTH»SINP/8VJ 1005 

COSTHsSQHTC1«»SINTH**21 1000 

1005 

SOTsO.O 1006 

CDTSO.O 1007 

I3al 1008 

CVMitO, 1009 

SVCneCVCllaO. 1010 

1011 

DO 1 NBg,J 1012 

1015 

C VELOCITY COMPONENT IN DIRECTION OF FLOW AT EDGE OF BOUNDARY LAYER lOlU 

SV(N)H(UH 2 fN)i»iCDSTH*li/' 3 QRTCRATl 05 +VH 2 rN 5 *SINTH* 8 INP)/aVJ 1015 

C ni8PUACKME.NT thickness IN DIRECTION OF FLOW AT EDGE OF B.L, 1016 

SDTaSnT+,5*(SVi:N)+3V(N-in*(ZCOMP(N)-ZCnMP(N»n5 1017 

r CROSS-FLOW VELOCITY COMPONENT 1018 

CV(N)BlJ*3INTH*(VM2(N)-UM2(N))/80HTCHATin) 1019 

r CROSS-FLOW DISPLACEMENT THICKNESS 1020 

CDTaCDT + . 5 *CCV(N)+cV(N-mAaC 0 MPINl-?C 0 MF(N-in 1021 

I8««IS 1022 

1025 

IF(A9S(CV(N)).cVM 51#1.2 102tl 

2 CVMaABSfCVfN) ) 1025 

1 continue 1026 

1027 

CDTaARSCCnT) 1028 

anT«ZCOHP(J).»8DT 1029 

1030 

RETURN 1031 

end 1032 


subroutine DMNSIONCS.UfCOSP.DELTAl jTHETAl ,ILP,CFX2.CFY2,CDFX, 1035 

ICDFXINF) 103« 

1035 

C calculates DIMENSIONALISING factor and dimensional 8.L» thicknesses, 1036 

1037 

CnMMnM/TE8T/RNLCln),INTRL,IEP 1038 

1039 
1 0 U 0 

DO I Npl jINTRL lOAl 

loa? 

C iiCALE REYNOLDS NUN8ER TO STANDARD FORM, 10«5 

IFClFR.EO.nRNBRNL CM) 10a4 

IECIFR,EQ,2)RNaRNL(N)/C08P lOiiS 

IF(IFR.EQ,3)RNaRNL(N)/C0SPw*2 10A6 

0a80RTCS/CU»RN*C0SP)) 10^7 

D1«0*DELTA1 10A8 

02PD*THETAl 1099 

IPCILP.Ea,0)G0 TO U 1050 

WRITE(2,S)RNL(N) ,n,Dl ,02 lOSl 

3 FORMATCIX, ITHREYNOLDS NUMBER* , F 1 0 , 0 » 2X , 23H (DTMEN8 lONAL Z)/CHORO«) 1052 

l,Fe.6,2H)Z,2X,10HDELTAl/C« ,F8,6,2X,10HTHETA1/Cf ,F«,6) 1055 




GO TO 1 10S« 

a K«lTP(2,5)RWL(Nf5,0,ril,D2,CrxH,crY2,C0Fy,Cr!rXINF iC55 

5 FOHMAK IX, ITHREYWOLOS KiUMBER* , F 1 0 , 0 » 2X , 23 M(DIm£('J 8 lONAL Z)/CHORO«) 1056 

l,FR.fc,2H)Z,2X,lCHDELTM/C«« #FR,6, 2X , J OHTHET 41 /C6 .FB.fc/ 105T 

25H CFX«C10,5,?xaHCFY*EiO,3,2X5HCOfXBE1 (l'.5,2XeHCOFXINFii€10,33 1058 

I CONTINUE 1059 

1060 

RETURN 1061 

ENO 1062 


1065 

SUBROUTINE GEOMTPY (l,RHO,CH) I0b« 

106? 

C given aerofoil CO-OBDINiTES X AND Z, TRANSFORMS X TO THETA, COHPUTES 1066 

r DISTANCFS around SURFACE SfX), NOSE RADIUS AND SECOND DERIVATIVES OF 1067 

C S(TMETA) FOR USE IN CUBIC SPLINf INTERPOLATIONS, 1068 

1069 

CnMMClN/r,E0M/XA(56?),ZA(56?) 1070 

C0MM0N/SFX/STHC565) ,TH(56S) ,F5TH(56?),INT<J,FZTH(36S3 1071 

DlMENSrCN SDU70) ,THED( 170) 1072 

1073 

Fl«3.lttlS9265 107<l 

INT«Mi»lNT4<il 107? 

1076 

C IS AEROFniL CAMBERED 1077 

IFCT.EG.O) GO TO 102 1078 

1079 

r SET UP LOteER SURFACE Cn-ORDINATES FDR SYMMETRICAL AEROFOIL 1080 

DO 103 N«l , lNT«Hj 1081 

XA(NJi»-XAC2*IWTO-nj 1082 

ZA CN)it*ZA (2*INTB-N) 10B3 

103 CONTINUE * 1080 

1085 

C TRANSFORM X TO THETA 1066 

102 NLF(lNTo-n*CItn 1087 

DC 105 -N«2,NL 1088 

TH(M)»THXCXACN5/CH) 1069 

105 continue 1090 

TH(1)«0, 1091 

1F(I,EO,0)CO TO J 1092 

1093 

THCIMTRJfPT 1099 

INT9b(1+1)*INT4«I 1095 

1NT9M16IMT9-1 1096 

1097 

I THCINT9)62*PI 1098 

1099 

DO 105 N«l, INTO 1100 

2A(N)bZA(N)/CH 1101 

1 05 CONTINUE 1102 

1103 

C COMPUTE INTERPOLATING FUNCTION FOR ZCTHCTA) UOO 

CALL C8S(TH,ZA,FZTH, INTO, 0,0, 0,0) 110? 

1106 

tF(RHO,GT.O, ) GO TO 108 1107 

1106 

C COMPUTE NOSE RADIUS IF NOT SPECIFIED 1109 

CALL C8I(TH,ZA,FZTH,INT0,PI,R0T.RH0) illO 


95 





1 n n 


RHas2*(WH.0i**2) nil 

1 112 

1P« CONTINUE 1115 

COMPUTE A«C LENGTH STH(N) TO EACH XA(N),2A(NJ AIRFOTU POINT, llU 

STHCN) is MFASURFn FROM LOWER SURFACE TRAILING EOGE TO UPPER 1115 

SURFACE trailing EDGE ANO IS APPROXIMATED AS THE CHORDAL 1116 

r DISTANCE PETViEEN aTHFOIL COORDINATES. 1117 

c ina 

1119 

STHdiaP, 1120 

DO 111 N*a,INTA 1121 

m STH(N))“STHCN-n+SDRTC t(XA(N)-XAfN-n )/CH)**2+tZA(N)-ZA(N-l nw*23 1122 

1125 

C COMPUTE interpolating FUNCTIONS FOR 6CTHETA) AT THETA(X) 112R 

CALL CSGCTH,8TH,FSTHnNTU,O.O,0.Ci) 1125 

1126 

RETURN 1127 

end 1128 


1129 

SUBROUTINE IFPRUTnFPT,lLP,LA3T,JACKPOT) 1130 

1131 

C OETFRMIMFS PRINT out REQUIRED 1152 

C BOTH PARTS-JACKPOTBli FIRST PART«UACKPOT»5, SECOND PART- J ACKP0T»2 1155 

C NO PRINT DUT-JACKPOTbO 115« 

1155 

c HAS non-COnVFRGEMcE nccUPHED 1136 

IFCLA8T.ED.23GO TO 2 1137 

1 1 36 

C has LAST POINT requested BEEN COMPUTED 1159 

IFCLAST.EQ.nJACKPOTiel lUO 

1191 

C IS FULL OUTPUT REQUIRED AT EVERV POINT OR AT THIS PARTICULAR POINT 1192 

'IFCIFPTUo.I.PR, nFPT.Se,2.AND.ILP,EQ'.mjACKPOT»l 1193 

1199 

C 15 FULL OUTPUT REQUIRED AT LISTED POINTS AND THIS IS NOT ONE Of THEM 1195 
IFCIFPT.EQcE. AND.lLP.EO.O)JACKpnT«3 I19fe 

RETURN 1197 

1198 

C HAS SEC&nD PART OF PRINT OUT skipped AT END OF LAST SUCCESSFUL STEP 1199 

2 lF(JAeKP0T,ED.5J60 TO 3 1150 

JACKPOTpO 1151 

RETURN 1152 

1153 

5 JACHPOT*2 1159 

RETURN 1155 

END 1156 


1157 

SUBROUTINE INSTABtC08P»8,U»CH,RATlQ, AMESO/AWINFSDJ : 1158 

1159 

C EVALUATES THE CROSS-FLOW REYNOLDS NUMBER.. CHl’. 1160 

1161 

COMHDN/TEST/RNL(10),INTRL>IFR 1U2 

CnMMC!N/CR088V/EVn70),CVn 70 3 ,SOT,CDT,CVM U65 

1169 

he - 


r CALCULATE CHI/SOPT (REVnOLDS NUMBER) 

C CROSS FLOW REYNOLDS NO, T» BASED ON THE fiTMTMUM KINEMATIC VISCOSITY 

r coEFFiriEMT • either free stream or fore value 
FACTM ,/RATI0F*1 ,5 
IFCA^’E3D,5T, AMINesd) FACT«1, 

CnRR»COT*Sf3RT(CnsP*5*CH/tl)*EACT 

WRITE(a,3) 

5 FDRMiTnx.SOH*** SWEEP INSTABILITY TEST ***) 

DO 1 N«1,INTRL 

C SCALE REYNOLDS NUMBER TO STANDARD EDRM, 

IFtIFR,EQ,nPN*RNL(N) 

IFCIER.EO^gJRNsRNLCNT /CDSP 
IF tIFR,E0.5)RNBRNL(NT/C08PAw2 

CHl*CDRR*SnRT (RN) 

WRITE (2.4) RNL(M) , CHI 

4 FORMATIJH ,17HREVNoLDS NUMBER* , F 1 1 , 0 , 2X , 1 RHCHJ (OWEN-R AND ALL) ■ *F7 
1 , 2 ) 

i CONTINUE 

RETURN 

end 


1165 

1166 

1167 

1168 

1169 

1170 

1171 

1172 

1173 
1179 

1175 

1176 

1177 

1178 

1179 

1180 
1181 
1162 
1183 
1169 

ties 

1166 

1 1 67 

1168 


SUBROUTTWE PLI S T ( IFPT , NL I ST , SX V , SX VINC , L. D8 , INT3 , CH , X ATT , DX , INC ) 1169 

1190 

C PREPARE LIST OF POINTS WHERE FULL OUTPUT IS REOUIREO 1191 

1192 

COMMON/OPLIST/OPX(2O0) iQP8(2P0) 1193 

dimension 8XVn),DUHp(365).Xt2),8(2),8XVINC(l) 1194 

1195 

IF(,TFPT-5)301 , 303.309 1196 

301 CALL STHFRMX (NLIST.OPX, OPS, DUMP, D, INTS.CH.XATT.H) 1197 

DO 310 N*a,ls!LIST 1198 

OPS(N-l,)«OPS(N) 1199 

310 CONTINUE 1200 

GO TO 305 1201 

1202 

303 NLI8T*L*1 1203 

DO 306 N«1 ,NLI8T 1209 

OPS(N)b8XV1NC(N+1 ) 1205 

306 CONTINUE 1206 

GO TO 307 1207 

1208 

509 N«l 1209 

IF(INC,EQ,0) GO TO 308 1210 

309 X(2)*N*0X 1211 

CALL STHFRMX ( 2 , X , 8 , DUMP , D. I NTS, CH , X ATT . 8) 1212 

0P8(N)hS(2) 1213 

IF(0PS(N).GT,8XV(L))G0 TO 307 1219 

nlISTpN 1215 

NpN+i 1216 

GO TO 309 1217 

308 DO 511 N«l,L 1218 

OP8fN)BFLOAT{N)*DX 1219 

IF(0FS(N) ,GT,8XVTNCCL)) go to 307 1220 

511 NLISTpN 1221 


47 



1222 

1223 

122tt 

1225 


1226 

1227 

8UBffDUTINE PR! NT (X , 8 ,U x DU , J , DZ , INT 1 , J ACKPOT » PBI i LC , ANGLE2) 1228 

122P 

C0N^'0N/RESIiL7SXl»’H2(i7f') ,LlN2t 1701 ,VM2 Cl 7n) , CELT A1 X THETA 1 xNOxOUOZ 1230 

IxPVDZ 1231 

COHMON/CR088V/'SV(1701 ,CV(1701 xBOTiCnTxCVM 1232 

COmmON^COMPRES/INC, A«1NE3Dx AMESD,SINP,C06PxGAMHAx 6AMI ,GAh2xCAH3x 1233 

IGAN«i, AMFS3DiRATIOxSCDHP,ZCOHP( 170 3 xT(170 J ,RH0D{170) 123ti 

niMFMStON XPRINTClol , 1235 

DATA XPPlNT/,l,.2,.3x‘.ttx,Sx.6x,7x.8,,P,'.9P/ 1236 

1FCJACKPOT.EQ.21 GO TO 7 1237 

1236 

v^RITECPxlll X,S,SCOMP,U, AMESDxOUxNO 1239 

It FORNAT(lHn,3HXK , F9 . 6, i)X , 3H8« , F9 , 6 x RX , 7H8COHP. » F9 . 6 . «X , 3HU« , 12R0 

1P9*6,UX7HAHE3D« x PR . 6 x ilX 1 1 HOU/D ( 8/L 1 P , PI 2 , 6, ilX 12, I IH ITERATIONS) 12«1 

12R2 

ANG«ieo./3,l«159g65>»ANGLE2 12ttS 

]P(LC.NE,01WRITEt2,2O) 1244 

IFaC.EQ.2)i«RlTEC2i2Ci) 1245 

20 FClRNAT(tHO,40HSTEP-LENGTH HALTED AFTER NriN-CONVERGENCE 1 1246 

wPITEf2,12)DELTAl ,THETA1 ,PUDZ,DVDZ,ANP, 1247 

jg FORMAT! IXxftHDELTAlP , FP , 4 , 4X , 8HTHET A 1 • ,F6, 4 , 4X # 1 2H (OU/OZ 1 Z *0 a, 1248 

lF6.4,2Xl2HCDV/D2)ZaO bE 1 0 , S , 2Xl 4HAIRF0IL 8L0PEaEl0,31 1249 

IFCJACKPOT.EO.Sl GO TO 8 1250 

C 1251 

r. TEMPORARY FIX ON PRINT OF VELOCITY PROFILES AT EVERY lOPC CHORD 1252 

C 1253 

IFCS,LT,l.E-oa) GO TO 7 1254 

C 1255 

rFCX,Lt.,n05) KbI 1256 

IFCX,GT,XPRlNT(i<n GO TO 7 1257 

GO TO 8 1258 

1259 

T IFCP6I,LT..0001) GO TO IS 1260 

V.RITE(2xt5)80T 1261 

15 F0RKAT(1H0x36H8TREAM flow displacement THlCKNESSa xPlO, 6) 1262 

WRITE(2x 161CDT - 1263 

16 FORMATCIHO^SSHCROSS FLOW DISPLACEMENT THICKNESSa ,Ft0',6) 1264 

WRITE(2,i9)CVM 1265 

19 FORMAT(1H0x26HMAXx CROSS-FLOW VELOCITYa xF10,6) 1266 

13 WRTTEC2, 14) . 1267 

14 FORMAT(4XlHZ,lOXSHZCOMPrexiHW,tlXlHU,11XlHV,lOXSH8TV,lOX3HCFVi 1268 

19XlHTx9XaHRH00) 1269 

1270 

17 NlalNTl+1 1271 

DO 2 N«l, J 1272 

IFCN,E0Vl«OR,N.EO.2) GO TO 3 1275 

IFaNTl.GT.2.AND,N,LE.lNTl)G0 TO 3 1274 

IFCN.EO.Nl .OR’.N.EQ. J)GO TO 3 1275 

GO TO 2 1276 

1277 

3 ZaN*DZ-DZ 1278 


307 IFPTag 
505 RETURN 

end 


1*8 


IFCP8I.lt. .0001)60 TO 1 1279 

W«ITE(2»<>)Z»ZC0MP{W) ,WM2(M) ,UM2(N) , VH?(N) ,8V(N) ,CV(N) ,TCK‘) *«HOD(M) 1280 
GO TO 5 1281 

1 '«PITEC?»9)Z»ZCDMP(N),^M2(M),UM2(N).VM?Cn),TCN),RH0D(N) 1282 

6 FORMAT(9CF10.8,2y)) 1281 

9 F0RMAT(5{Fl0.fc,2X),2-4X,?CF10.6,2X) ) l2Btt 

5 IFCN.EQ.NDMieMl + lMTl 1285 

2 CONTINUE 1286 

IF (S.LT, I .E-OB) GO TO 8 1287 

K«K+1 1288 

1289 

8 RETURN 1290 

END 1291 


1292 

SURROliTINf RELAMCU,8.&INp,C08P,KFAX) 1295 

129U 

C evaluates the RE-LAM1NaR!8aTTON PARAMETER. K, 1295 

r CHAS NOT been modified for COMPRfSSlBTLTTV) 1296 

C 1297 

1298 

dimension UC«),S(2).A(2) 1299 

COHMnN/TEST/RNL(in),IMTRL.lFR IlOO 

REAL KMAXnO) 1101 

1502 

WHITFC2,u) 1501 

<1 FaRMATdX.SOHAt** RELAMINaRISATION check ***) 130<» 

1505 

DO 1 Nssi ,2 1306 

A(N)«8QRT(UCN)**2+8INP**25 1507 

1 CONTINUE 1508 

1309 

08pS(2)*S(n 1510 

DlJ«U(2)+Utn 1311 

RCDSPK»0*CAC2)-An))AOU/(O8*(AC2)+A(l))**3) I5l2 

1515 

DO 2 N»1 .INTRL 1.319 

CAYaRC08PK/CRMLCN)*C08P) 1515 

TF(IFR.E0.2)CAYPCAY*CO8P 1516 

IFCIFR,Ea.3)CAY«CAY*COSP*K2 1317 

1518 

IFtCAY.6T.KMAX(N))KMAXCN)«CAY 1519 

1320 

HRITE(2.5)RMtCN),CAY,KMAXCN) 1521 

3 FORMATCIH ,17HREYN0LD8 NUMBERW ,F10.0,5X,5HK» .E10.3,3X,6HKMAX» ,E 1522 

110.31 1525 

2 CONTINUE 1529 

1525 

C IF K.LT.1/2 K(MAX) set parameter TO AVOID COMPUTING KOOMNSTREAM. 1326 

IF CKMAxaNTRL).GT.2*CAY)KHAXnNTRL)»-1.0 1527 

1328 

RETURN 1329 

END 1530 


1331 

SUBROUTINE 8TPLNTH tSNEXT , 8, INTHOLD . D8. D81 » D8Z, NgXT . Nl.1 8T.I AST , 1532 

49 



lIFPT,LC,ILP»U,U8TFP,0llDS,X,SH,ITC,ANGUE2,ViF) 133X 

133U 

f CALCULATES LFMGTH OF ME^T STEP. 1335 

133E 

CnMMON/XSANDU/UMfS*)5 5 ,THXVC3655 , FUTH ( 3(,«! ) , XV (3fc5 ) , CPUK ( 365 ) » 1337 

1SXVC365) »3VVI)gC(365) ,FSVSINC{365) .L,SATT, INT3#CH,I8P 133B 

COMMON/nPLlST/OPX(?00^ .OPSC200) 1339 

CO^‘Pn^4/cOPPRF8/I^C,AMI^F3D, AME3D,8IMP4CnSP,GAMHA ,GAM1 ,GAM2,GAM3» 13«0 
lGA‘^U,AHFS3PiPATin,SCOMP,ZCnMPtl70),Tn7O).RHOO(17Q5 13«1 

OP^FnSION UC<n,S(2) 13«E 

13«3 

ILPpO 13«« 

LlNFARaO 13«5 

I3«b 

C A'AS LAST STEP SUCCESSFUL 13«7 

IF fUC,Ew,0)GO TO 9 13i»6 

I3i<9 

C WAIVE STEP UEMT,’'TH AFTER N0N»C0NVER5EA!CE 1350 

1 IF (081 ,GT. aOlikDSlGO TO 1 1 1351 

f STEP LEK'GTH less than MINIMUM PERMITTED- END CALCULATION, 1352 

LASTb? 1353 

GP TO 8 I35^i 

1355 

11 £)S1b081/2,0 1356 

LASTbO 1357 

GO TO 7 * 1358 

1359 

9 IF(SNEXT,LT,0.) go to ? 1360 

S(23aS»FXT 1361 

D8lBSC2)-8(n 1362 

SNEXT«»1,0 1363 

GO TO « 136a 

1365 

2 IFdNTHOLO.EO.OlGO TO 35 1566 

r value of next point has been held while a listed output point was 1367 

r COMPUTED, 1368 

D3U0SH 1569 

C step length- is held at D8H FOP INTHOLO STEPS, 1370 

INTHOl 0«tNTH0L0-l 1371 

GO TO 7 1572 

1373 

c standard step length, i3?a 

33 DS1«D8 1375 

C CHECK ratio (proposed LENGTH DF NEXT 8TEPI/(LENGTH OF LAST STEP) 1376 

« IFCDSI ,LT.1.2*0S2)G0 TO 7 1377 

1378 

C KEEP STEPLEMGTH CONSTANT OVER A NUMBER Of 8TEP8 DEPENDING ON ABOVE, 1379 

INTH0LDarNT<5*DSi/D8Zl-5 1580 

IF(INTH0LD,GT,53 INTHOLDaS 1381 

1382 

C limit length of NEXT step TO TWICE LAST STEP. 1363 

IF(lNTH0LD,E0.5)08r!S2*DSZ 1384 

OSHbDSI 1365 

INTHOLDBINTHOLP-I 1386 

* 1387 

C PROPOSED VALUE AT ENO OF NEXT STEP. 1388 

7 S(25«S(1 ITDSI 1389 

iFdFPT.EO.n GO TO 10 1390 

C IS PROPOSED VALUE OF X LESS THAN THAT of NEXT LISTED OIITPIIT POINT ll91 

IF(S(2).LT,0PSCNE)(T)) go TO to 1392 
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o o n 


15«3 

r WCPLACE proposed value with that CIE NEXT LISTED OUTPUT POINT 139« 

S(2)«OP8(NEXT) IS95 

I39fe 

C IS this the last POINT IN output LIST 139T 

IE(NEXTtI.EO.NLIST)LAST«1 139S 

1399 

C ADVANCE OUTPUT LIST COUNT l<JOO 

next»next+i laoi 

1LP»1 laoz 

IFaNTHOLO.NE.OISNEXTBSdTADSH ia03 

laoii 

C IS PROPOSED VALUE OP X LESS THAN THAT OP LAST VELOCITY DATA POINT l<f05 

10 TP (Sf2) ,LT, SXVINC(L)-. 0001*01) GO TO «a 1«06 

laoT 

C VALUES AT LAST VELOCITY DATA POINT, 1«08 

S(?)«SXVINC(L) 1R09 

’^■XV(L) 1«10 

u(2)«u“.tL) i«n 

()UDS«CUMCU-UM(L-1 ))/(8XVINC(L).8XVINCCL-1 n 1^12 

LA6T»1 lUll 

60 TO 45 1414 

1«15 

U« IP (1NC,E»,0) 60 TO 4#) 14ifc 

SC0HP«S(2) 1417 

GO TO 47 1418 

C 1419 

C INTERPOLATE SXV(SXVlNc) TO PINO LOCATION |N PHYSICAL PLANE, 8C0HP, 1420 

CORRESPUMDING TO LOCATION IN TRANSPORMffO, INCOMPRESSIBLE PLANE, 1421 
S(2) 1422 

1425 

4fc CALL CSI CSXV1NC,SXV,FSVSINC,L»8(2) ,8COMP,ROT) 1424 

C 1425 

C find velocity AT END OF PROPOSED STEP 1428 

U7 CALL XNnpRM8(SC0MP,UC2)»rJUDS,X,lTC,ILP,ANr.LE2,LlNEAR»THETAS) 1427 

1428 

r nTD ncim-cOnVerrENCE OCCUR IN SUB-ROUTINE XNDpRMS 1429 

IP(ITC,Et3, 20)60 TO 15 1430 

l«31 

nSl»8C2)-8(l) 1432 

1433 

C check THAT USTEP IS NOT EXCEEDED AND REDUCE BTEPLENSTH IP NECESSARY, 1434 

1435 

17 IP{AB8fU(2)-Um).LT,l)STEP)G0 TO 19 1436 

IP(ILP,EO,l)NEXV«NEXT-l 1437 

ILP«0 1438 

1419 

C ITERATION TO FIND 8 FOR CU(l)+U8TEP) 1440 

INTLIiO 1441 

St2)P8(?)-DSU(i-U8TEP/(U(2)-U(l n) 1442 

IFtINC,EO,0) GO TO 48 1443 

SC0HP«8(2) 1444 

GO TO 12 1445 

48 CALL C8nSXVlNC»SXV*F8V8INC,L*8(2),8CnMP,R0T) 1446 

12 CALL XNDFRn8(9C0MP, 11(2), duds, X,ITCiILP.AN8LE2, linear, THETAS) 1447 

IF(ITC,E«,20)GO TO 15 1448 

IF(A88(ll{2)-UCl)-U8TEP).LT.,01*U8TEP)r,0 TO 111 1449 

INTU«1NTU+1 1450 

IF(INTU,E»,25) GO TO 24 1451 

S(2)R8(2)-CU(2)-U(n«U8TEP)/DU08 1452 



IFU^C.EO.O) GO TO ue l«53 

SCOMp«S(E) HSa 

GO TO 12 

1^56 

15 i“WITC(2,26) 1«57 

2fc FDSMAT(lHO,«lHMON-COKiVePGENCE l^J X *NO tl FROM S ROUTINE) 1«56 

GO TO 120 l«59 

1U60 

24 IF(4BSCLi(2)-Ufn).GT,U3TEP) wRlTE(2#3) 1461 

3 FOP’^ATCIHO,* INCREMENT In L) HAS EXCEEREO SET LIMIT*) l«62 

1463 

r either non-convergence has occurred in finding thetas for a given 1464 

C SCOHR in S/R XMDFRMS or non-convergence has occurred in finding U(2) 1465 

f that satisfies TME USTEP criterion, set ll(2)«lin) ♦ USTEP and use 1466 

C linear interpolation to find S(2), 1467 

1466 

120 LINEARbi 1469 

U(2)«Um+USTEP 1470 

DO 60 N«1,U 1471 

IF(U{2),LT,UM(NJ)G0 to 65 1472 

IF(8XVINC{N).GT,.5)60 TO 63 1473 

60 CONTINUE 1474 

63 PRINT 66 1475 

66 FORMAT(/H- LEADING-FDGE LINEAR INTERPOLATION OF U(2) VS, 8(2) TO 1076 

IFIND seal FOR A RIVEN U(2) RESULTS IN 8(2) GREATER THAN ,5*) 1477 

STOP 6666 1476 

65 FACTb (U(2)-UH(N-l))/(UMCN)-UM(N-n ) 1479 

S(2)avSXVIMC (N-n+ (SXVTnC (N)-8XVINC(N-l ))*FACT 1460 

THETAS*THXV{N-n + (THXV(N)-THXV(N-n)*FACT 1481 

DUDSs(U(2)-Uf 1 n/f8(2)-S(n) 1482 

IF(INC,rC3,n)GO TO tin 1483 

SCOMPb6(2) 1484 

GO TO 112 1485 

rm CALL CSI(8XVINC,SXV,F8V8INC,L»8(2),8C0MP,RDT) I486 

112 CALL XMDFPM8(5C0MP,U(?) *DUD8,X,ITC, ILPr*NGLE2»LlNEAR, THETAS) 1487 

111 LASTaO 1488 

1P89 

r LENGTH OF NEXT STEP, 1490 

19 DSlsS(2)«*8Cn 1491 

r S LOCATION OF next STEP 1492 

i?H»HF*S(2) + (l.-wF)*8(l) 1495 

1494 

e return 1495 

END 1496 


1497 

SUBROUTINE STHFRMX(J,X,8jTHXV»pS07,INT5.CHiXA7Tr8ATT) 1498 

1499 

r FINDS 5(N) AT POINTS X(N) FOR W*1(1)J AND D8/DTHETA AT X(l) 1500 

C where 8 18 measured FROM THE ATTACHMENT LTNf, 1501 

1502 

1503 

COMM0N/8FX/8TH(365),TH(365),F8TH(S6S),lNTfl»FZTH(365) 1504 

dimension X(l), 8(1), THXV(l) 1505 

1506 

XCDbXATT 1507 

DO 1 NbI , J 1508 

THXV(N)bTHX(X(N)7CH) 1509 
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IF(lKiT3.e0.05Ga TO 2 I5l0 

1511 

CALL CSI(TH,5TM,F8TH,int«»THXVCN),8(M),DI0TH) 1512 

1513 

TF(M.GT, DGO TO 3 15lt 

OSDTbDSDTH 1515 

GO TO I I51fc 

1517 

g S(M)«X(N)/CH 151ft 

GO TO 1 1519 

1520 

3 S(N)»S(M}-S(U 1521 

1522 

1 CnwTIHJF 1523 

SATT»S(n 152U 

8(n«0, 1525 

1526 

RFTLiftNi 1527 

FnD 152ft 


1529 

ftllRSnUTIMF TRA^JS(8C»USSM,THETA1 ,llTWn,LIONE, 1ST, JACKPOT) 1530 

1531 

r F8T1^1ATE8 THE POSITIONS OF VlftCOUS TMSTABTLITY AND SUBSEQUENT 1532 

C TRANSITTOM, 1533 

1534 

COMNfOKi/TEST/RNLflO),TNTRL.IFR 1535 

COMmp^/sUPTRAM /R0S(17)» AM(17)»GRAN(15). AMT(13) 1536 

rOHRON/COMPRES/U'C, ami NFSOjARESO, 81 NP,CnSP, GAMMA, GAM 1,GAM2,GAM3. 1537 
1 GAM(j,amFS 3D, RATIO, 3COMp,2COMP(l70),T(170),RHOnn70) 1538 

COMmoN/YSANDU/UM(365),THXV(365) ,FUTH(365),KV(S65) ,CPUMC365) , 1539 

ISXV(365),SXVINCC365),F8V8INC(365),L,8ATT,TNT3,CH,I8P 1540 

0IMFV8I0H SCTW( 10),SCI tl 0 ) , 8UMM ( 1 0 ) , RTCL ( 1 0 ) , POL Cl 0 ) , 9TL ( 1 0 ) , 15«1 

IRTHIO) ,sciCOMPCin) I5a2 

1543 

B*(l,/RATIO) a*GAM1*S0RT( 1,+AMF830aCQ8Pa*2)/(1 .♦AMF83D*UTW0a* 2)**1 , f,3aU 
15 1545 

IFCIST.NE.l) GO TO 2 I54i 

1547 

C 8FT certain stores TO ZERO WHEN 8/R IS ENTERED FOR FIRST TIME, 1548 

8CL»0,0 1549 

SCOMPLbO.O 1550 

RATIOLbI ,/( l .+AHFR30*C08P«*2) 1551 

00 4 N«1,INTRL 1552 

8CTR(N)»-0.0 1553 

8CIfN)«*0.0 1554 

8CICOMP(N)«-0.0 1555 

8UMM(N)bO.O 155* 

a continue 1557 

155ft 

ISTbO 1559 

1560 

C evaluate CLAM0OA)2 - 8A8E0 DM MINIMUM KINEMATIC VISCOSITY COEFFICIENT, 1561 
C ■ EITHER FREE STREAM OR EOGE VALUE 1562 

FACTsl ,/RATIOA*l .5 1563 

IF(AME3D,GT, AMINF30) FACT»1, 1564 

2 EM«9*ftC*U8ftM*THETAl A a2/UTWQ*FACT 1565 

kiRITE(2, 10) 1566 
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in FORM-1, TUX, 35H*** TRAMSITION TEST (GRANVILLE) ***) 1567 

DO 3 N«1,INTRL 1568 

WR1TEC2, inRNL(N) 1569 

11 FORMATClH ,17HREVMnL08 NUHBEH* ,Ft0,01 1570 

IK(SCTR(N) .GT’.o'.) go TO 1571 

1572 

C SCALE REVMOLDS NUMBER TO STANDARD FORM, 1573 

IFdFR.ED, URNsRNLCN) 157<» 

IF(IFR.EQ.2)PN*RML(N)/C0SP 1575 

IF(IFR.EG,3)RNiiRNLCN)yC0SP*Aa 1576 

1577 

C EVALUATE R2 1578 

C REYNOLDS NliMRER IS RASED ON mimmuh KINEMATIC VISCOSITY COEFFICIENT* 1579 
r, either free stream or edge value 1580 

RD*80RT(RNi»SC*UTwniirC08P/RATlO)*THETAl*FACT 1581 

1582 

wRITE(2v22) RD 1583 

22 FDRMAT(1H+,30X,7HRTHETA*F6.1) 15841 

IFCSCTRtN) ,CT.O,) GO TO ]U 1585 

C HAS INSTARILITY REEN PREDICTED UPSTREAM OF THIS POINT 1586 

TF(SCICN).GT.0,1 GO TO UO 1587 

1588 

r FIND CRITICAL VAI.UE OF R2 FROM STUARTS CURVE, 1589 

DO 6 J»2, 17 1590 

IF(AM(J),GT,EM)G0 TO 7 1591 

6 CONTINUE 1592 

J»17 1593 

I594i 

7 AB(RnS(J-n-KR08( J)*R08( J«n )*(EM»AM{,T-n)/(AM( J)-aM( J.l ))) 1595 

RTC«in,iir*A 1596 

IFCRO.GT.RTCIGO TO 9 1597 

1598 

RTCUN)*RTC 1599 

RDL(N)bRD 1600 

IF{ JACKPOT, EP,3)G0 TO 3 1601 

WRITEC2, 1101RTC,EM 1602 

lin FnRMATClH+,4iiiy,UH RTHETCRIT«F6’.1,2X,6H L AM2«F6 , 3 » 2X , 1 OHNO INSTABI 1603 
ILITY) 160U 

GO TO 5 1605 

1606 

r interpolate for values AT POINT OF INSTABILITY. 1607 

9 SCl(N)o8CL+(8C-8CU*CRTCL{N)-ROL(N))/(RD»RTC-RDL(N) + RTCL(Nn 1608 

IFdNC.EQ.O) GO TO 20 1609 

SClCOHp(M5«SCI(N) 1610 

60 TO 21 1611 

20 CALL C81(SKVINC,8XV,F8V81NC*U»8CI(N),8CIC0MP(N),RQT) 1612 

21 RTICN)aRTCLCN)+(RTC-R7CL(K))*(8Cl(N)-8CL)/{8C-8CL) 1613 

EMI«EHL+(EM-EML)*(SCIfN)«8CL)/(8C«8CL) 1614J 

SUMM(N)*0.S*(EM+EMI)*(8C-SCI(N)) 1615 

GO TO 8 1616 

1617 

lao DIMB0.5*{EM*RAnnAt*GAMl+EMLFRATTOL**OAMl)*(8C-8CL) 1618 

8UMM(N)b8UMM(N)+DIM 1619 

1620 

C EVALUATE (R2)T-CR23I 1621 

B RTMRTIbRD-RTICN) 1622 

1623 

C evaluate (LAMBDA)2 BAR 1629 

AMB«8UMM(N)/(8COMP«aclCOMP(N)) 1625 

1626 
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C FIND CRITICAL VALUE OF (R2)T-(R2H FROM SRANVTLLCS CURVE, 

DO 15 K»2, 13 

IF(AMTCK) .GT.AMB) RO TO 16 

15 CONTINUE 
K«13 

16 RTCMBTIbGRANCK-1 ) A ( AMR- AMT ( K- 1 VI «K6R AN ( K ) *CR AN (K- 1) ) / 

1 CAMT(K)-AHT(K-in 

l.a ■<R1TE(2, 13)SCICDMP(N] 

15 FORMAT(lH+,««y,20HIN3TABlLITY 'AT 8/C« ,F6.<0 
IF(SCTR(n),GT. 0.)SO TO 5 
IFCRTMRTIgGT.RTCMRtn GO TO IT 
WRITE (2, 19) 

19 FDRMAT(lH*,72y,l5HNO TRANSITION) 

RTL(N)«RTCMRTI-RTMRTr 
GO TO JOO 

17 SrTftCN)HSCOMPL+RTL(N)*C8COMP-8COMPL)/(RtMRTI-RTCHRTI+RTL(N)) 

5 wr1tE(2,18)SCTR(N) 

IR F0KMAT(1H*,72X,19HTRAN8TTT0N AT 8/C* ,F6,<1) 

100 WRlTF(2,10n RTKN) ,«TCMRTI,AMB 

V<H--E^Q-SjiUXi.r7M IN8TAP. RE’. Nn',*F6, 1 ,2K,8MRTC*RTI«F6, 1 ,2X,8 hLaM2BAR*F6, 


15) 

3 continue 


1627 

1628 

1629 

1630 

1631 

1632 
1631 
163R 

1635 

1636 

1637 

1638 

1639 

1640 

1641 

1642 

1643 

1644 

1645 

1646 

1647 

1648 

1649 

1650 

1651 

1652 


SCL»8G 

SCOMPLMCOMP 

EML«EM 

ratiol»ratio 

RETURN 

end 


1654 

1655 

1656 

1657 

1658 

1659 


1660 

BLOCK n.ATA 1661 

1662 

C STORE TABLES DERIVED FROM STUARTS AND GRANVILLE8 CURVES FOR USE IN 1665 

C SUB-ROUTINE TRANS. 1664 

1665 

COMMON/SUBTRAN /Rn8(l7),AM(17),6RAN{l?),AMT(l3) 1666 

DATA R0S/l,392,l,464,l,57S,l,T,l,84,2,0l6,2,2H4,2'.45e,2.7l5,2,956, 1667 

1 3. 1 55, 5. 3 10, 5. 452, 3. 57, 3, 676, 3, 734, 3. 768/, AM/., 06, •. 05, -.04,., 03* 1668 

2-. 02, -.01, 0,0,, 01, ,02, ,03, .04, ,05, ,06, ,07, ,08, ,09, ,10/, BRAN/450,, 1669 

5460, ,4S0.,504,, 548., 61 0,,706,, 836, ,1000’ ,1195., 1440',,1720.,2o46,/, 1670 

«AMT/-,035,-,030,-.025,-,02,-,015,-.01,-.OC5,,0,,005,,010,.015*,02, 1671 

5,025/ 1672 

end 1673 


1674 

SUBROUTINE VEL0CT8 flNTV,ca8P5 1675 

1676 

C COMPUTES U FROM DATA 1677 

1678 

CnMM0N/X8AN0U/UM(365),THXV(S65),FUTH(365),XV(l65),CPUM(365), 1679 

iSXV(365),8XVlNGC565),FSVSlNC(365),L*8ATT,INT3,CH,I8P 1680 




iu. 


OF POOR OtJARTTY 


16B1 

UPluL+l H>82 

00 I N»8,LP1 lt>83 

IFtlNTV.ED.SJGO TO 3 U8« 

1685 

t INPUT VELOCITIES WERE NON.DIMENSIONALISED w'.R'.T'. FREE STREAM VELOCITY 1686 
C perpendicular to leading edge, 1687 

UM(N)*UM(N) *COSP 1688 

GO TO I 168R 

1690 

C COMPUTE U FROM PRESSURE COEFFICIEN'TS 1691 

S CPUMCN)sUMfN) 1692 

UM(N)»80RTtC08P**H-CPUM(N)) 1693 

1 CONTINUE 169<< 

1695 

CPUMCn«C08P**(2 1696 

RETURN 1697 

end 1698 


SUBROUTINE. VGRAnAT(ALPHA,RH0,81,83,VGRAD,YV) 

C ESTIMATES DU/D8 AT ATTACHMENT LINE AND POSITION OF ATTACHMENT LINE, 


INCtl 

' XgJJ I-PHA .LT’. 0 . 0 n Nc «. 
ALPHA7rNC*Al,rlTA~ -r -., 
A«,017a55*ALPHA 


XATT« CTAN(A)*n+S3n**2/( ClT8n*A2A(TAN(A5*(lT83n*A2> 
VGBADaCOSC A)*( l + sn*nTXATT5/(RHO + 2*XATT) 

XV«-XATT*INC 


RETURN 

END 


1699 

1700 

1701 

1702 

1703 

1709 

1705 

1706 

1707 
i7oe 
1?^ 

1710 

1711 

1712 

1713 

1714 


1715 

SUBROUTINE HNDf RM8CS » U , 0U08| X , 2TC # ILP , ANGLE i LINE AR, THETAS) 1716 

1717 

c finds X(8) FROM 8 BY ITERATIVE METHOD# HENCE U(8) AND OU/D8, 1718 

1719 

C0MM0N/SFX/8TH(365)#TH(365)#F8TH(369),INT4,FZTH(365) 1T20 

C0MM0N/X8AN0tJ/UMt365).THXV{365),FUTH(365),XV(S65),CPUM(365)» 1721 

1SXV(565) »8XV1NC(365) ,F8V8INC(S6B) *L»8A7T,INTSiCH,T8P 1722 

1725 

C0MH0N/GE0M/XA(S65),ZA(368) 172« 

COMMON/COMPRES/INC, AMINF30#AME30# 8INP,C08P, GAMMA #GAMl,GAM2,6AM3» 1725 
\GAH4,AMFS3D#RATIO,8COMP,ZCOMP(170), 7(170), RH0D{170) 1726 

ITC«0 1727 

tF(LINEAR,EQ,nGO TO 15 1728 

1729 

1F(1MT3)1#1.2 1730 

1 THETA8pTHX(8) 1731 

D8DTH«*0,5*8INCTHETA8) 1732 

GO TO 3 1713 

1734 
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r. THETACS) FOR ESTIMATED 8 I735i 

? TEST»8+XVtn/CH 173hj 

lFrTEST,GT,i,5 GO TO 7 1737! 

TMETAS*THX(TESn 1738 

Cn TO tt 1739 

7 THETAStfe, 1790 

1791 

r FTND S AND DS/n(THETA) AT S 1792 

9 CALL C87(TH,STH,ESTH,INT9, thetas, 8X1, 08DTH) 1793 

TEST1«St8ATT-SX1 1799 

TF(ABS(TESTn .LT.n.ftOOOn GO TO 3 1795 

I79fc 

r I^-PROVE estimate FDR 8 AND EVALUATE THETA(8), 1797 

THETABBTHETAS+TFSTl/nSOTH 1798 

ITC«ITC + 1 1799' 

IF (ITC,LT,20)GO TO 9 1759 

1751' 

WP1TE(3,5) 1752 

5 FrRMAT(90HaNON-CDHVERGENCE IN S TO THETA PROCEDURE) 1753 

RETURN 1759 

1755 

C FIND U and DU/D(THETA) AT THETAC8) 1756 

3 CALL CSI (THXV,UM,FUTH,L»THETA8,U,DUDTH) 1757 

15 AHE20*AMINF3D*U 1758 

AVE5D*{AMEaO**2* (1 , + AHE83 D)A(AMTNesd*81NP)#i*2)/(i .♦AMP830*C08P**2) 1759 

AMEsn^seRT (AME3D) I7b0 

RATTO*n ,+.5*(GAMMA-1, )*AHESD**2)/(1,+AMFS3D) 1761 

DSCOSIpflATIO**GAMl 1762 

IFCLINEAR.Ea. non TO 20 1763 

DliD3P0U0TH/DSnTHii.D8CD8l 1769 

-- .?aJCpXTH(TWETAS)*CH 1765 

‘ ' — : 1766 

IF (ILP.EQ.t) RETURN — 1767 

CALL C8HTH,ZA,FZTH,INT9,THETA8,R0T»DZ0TH) 1768 

DX0Th«ABS(,5*81N(THETA»)) 1769 

IFCDXDTH.QT, r.0E-05)en to 6 1770 

ANGLE»3. 191592659/2. 1771 

RETURN . 1772 

6 ASLOPFPDZDTH/oxnTH 1773 

ANGLEBATANC ASLOPE) 1779 

RETURN 1775 

END 1776 


1777 

subroutine XSCPPNT (INTV) 1778 

1779 

C PRINTS OUT TABLE OF VELOCITY DATA, 1780 

1781 

C0MM0N/y8AN0U/UHt365).THXV(56f),FUTH(165)iXV(S65) ,CPUM(565), 1782 

18XVf365)#8XVlNC(365)<pF8V8INCC365)»L,SATT,TNT3»CH,I8P 1783 

1789 

DIMENSION 8(365) 1765 

1786 

IF(IN73.EQ,0)G0 TO 6 1787 

DO 5 Npl,L 1788 

8CN)f8XVCN)*CH 1789 

5 CONTINUE 1790 

6 CONTINUE 1791 



iT92 

IFCI'^iTj’,E(5.0.AMD.lMTV',LE,2)WRITE(2f n fXV(N),UMtN),Nlil,L) H’S 

1 F0RMATnHft,«X,2HXV,8X, IHU/CIH , 2 CPF . « , 2X ) ) ) IT94 

1795 

IPClMT3,Ea,0,AND,INTV',E0,3)wRrTP(2,.lJ f X Vf N ) , CPUM (N ) , UM C N ) , N»1 ,L ) 1796 

2 F0HMATnHn,«X,2HXV,7X,2HCP»9X»lHU/(lH ,S(F«,af2Xn) 1797 

1798 

IF(INT5.EfM,AN!5.INTV'.l.E,*)‘^PnKC2,5) (XVCN)»8XVCN)#SXVINC(N), 1799 

lUM(M),THXV(N),FMTH(N),FSV8l»^C(M5,N«l|tl 1*01 

3 PCIR|^AT(1H0,7X2HXV, 14X3H8XV# 11 X6H8XVIMC, 13X1HU,11XUHTHXV, 12X4HPUTH. 1801 

1UX7WP8V8IMC/(1H ,7F.16,«n 1*02 

1803 

IPCIIvTS’.EQ. I .*’'^0.TMTV'.EC3,3)WRITP{2.«HXV(N) ,8Xv(N),8XVINC(M» 18 04 

irP*!W(M) *UM(N) ,MBl ,L) 1805 

« FnRKAT(lH0,6X2HX\/,14X5HSXV,l3X6HSXVIMrinX2HCP,lilXlHU/ 1806 

IMH ,5P16.8n 1*07 

1808 

RETURN 1*09 

end I8l0 



APPENDIX E 


SAMPLE CASE 

The sample case consists of the computation of the boundary layer on 
the upper surface of a wing swept at 35° with the airfoil section shown in 
figure 2 subject to the suction distribution given in figure 1. This airfoil 
which is nominally 13^ thick was designed specifically for LEG use by Pfenninger, 
Allison, and Bobbitt using the inverse method in reference 6 to design the 
airfoil and the analysis method in reference 7 to modify the lower surface. 

The sample case free stream Reynolds mmaber is 11 x 10^, based on the chord 
measured perpendicularly to the leading edge. The free stream Mach number 
is 0.885 which gives a Mach number normal to the wing leading edge of 0.725, 
the same as the design value. The suction distribution shown in figure 1 
maintains laminar flow over the entire wing surface according to the criterions 
which were previously discussed. It should be noted however that no attempt 
was made to optimize this suction distribution; hence, it is expected that 
these suction levels can be reduced thereby reducing the skin-friction drag. 

The input for the sample case is listed below. The program prints 
this input as well as some computed quantities and that information is also 
listed below. A sample of the output is then shown with both the print out 
at a typical boundary-layer station and the boundary-layer profiles given. 

Figures 3 gives the distributions for this sample case of the x and y 
skin- friction coefficients along the surfiace from the leading to tpailing edge. 
This sample case required a total of I 8 seconds and 7^0 K storage for 
execution on the CDC CYBER 175 computer. 
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INPUT FOR SMPLE CASE 


50 100.00001 ,05 .01 ,05 1. 

I 
7 

. 

.5 .5 

• 

.5 ,7 

. . B 

.5 .<? 

. 5 0 . 

0 ,885063 l,a 

suction. TRiNSITIOM ANAYL8IS OF YAwgp V.ING LAMINAH BOUNDARY LAVER 

n. 

G 161 ,0003 
-.lOOOOOOOE+01 0. 

-.RRR«9851E + O0 51 8U75»6E-0i) 


• ,<>R60103RE + Qn 

«.RR 5555 feOE 4 .oo 
-.RR2139/JOEFOO 
-.9677638«E+00 
-. 982 A 2976 E +00 
-.97613759E+00 
-.96688912E+00 
-,960fc8937E400 
-.9515U675E+00 
*.9«107283E+00 
-, 930 ti 85 n 2 E 400 
«,<5ifi59803E400 
-.ROSSilSSAEfOO 
-,8922aft«OE400 
",87785050E400 
-,8626B915E400 
-.89679335E+00 
-,S3019369E+00 
-.81298151E400 
-,7953t|209E400 
•,77750733E400 
-,759675il6E + 00 
-,7i»197277E400 
-.72a«2ft86E+O0 
•.70699660E+00 
-.68963086E400 
-.67229757E4Q0 
•,65097195E400 
-,6376P323E400 
•.82021565E+00 
•,60273999E400 
•,5851931ZE+00 
-,56756606E400 
•.5«98«185E400 
•,53200957E400 
-.51908S35E+00 
•.49610193E400 
-,a701OTl6E+OO 
-.^ 8013100 E 400 
-,a«219968E+00 
-,tt243395«Et00 
-,il0657596E400 
-,Sfi«9358ftE+00 
•,371tt0931E+00 


B 12012006E-0S 
.269963U7E-03 
.06l6295iiE»03 
,69008il30E«»03 
,9«fJ51567E-03 
,1211 1269E-02 
.1<^722551E»02 
.17070B01E-02 
. 18927079E-02 
.2noui255E»02 
■.20136525E-02 
.18915196E-02 
,1607a39A£-02 
;11289977E-02 
,ii23136305-03 
‘.538795aiE-03 
,179} 06031-02 
,3^066190E-02 
.59978826E-02 
.01457877E-O2 
■. 1 1330604E-01 
.1«9S2555E-01 
.1S770622E-01 
'.226518971-01 
,2605R1U«E-01 
',3005B77«E-01 
',33468083£-01 
,36620239E-01 
'.39a84l50E-0l 
J.42054396E-01 
.ttU3S3053E-01 
,a6316717E-01 
^« 799 <J 70 aE -01 
,«938854«E-0l 
,505ll«9BE-01 
'.5H09030E-01 
.5P109325E-01 
.52619291 E-01 
■.52937205E-01 
,53 ft 63 loss- 01 
,5P99B775E-01 
’.52746910E-01 
,523114571-01 
.51698406E-01 




IS 

Y 


6o 



,S5<H5«9feE400 
. 33709 Q 61 E +00 
.3203a595E+00 
,S03fl5«9fcE+00 
,?e7TOaU6E+00 
.27lB<1737E + 0n 
,2!»6E13l7E + 0f) 
.2ao7iesoE+on 

»2098799fcE<^00 
,l94a5512E+00 
. 1790553 ae +00 
. 16578959E+00 
. lilP790jPE+On 
,13U17989E+00 
.120069fcaE+00 
.10t>5527fcE+0n 
,93700565E-Ot 
.B1558093E-01 
.70180912E-01 
.?9fel7866E»0i 
,il9B9l09UE-0l 
,U1009359E-01 
,52978«23E-M 
.2S799979E-01 
,19a79093E-01 
.H01237fcE«01 
,939fil66CiE-02 
,56^E2706E-02 
.2S150762E-02 
.Qa968gil2E-03 
. 30 aB 596 UE- 0 U 

’,850302ii7E-03 
„l7782539E-02 
.3U83no7E-02 
,56B29827E-oa 
,B5fe3573fcE-02 
iU07l07?E-ni 
.182U803E-01 
.2O98B30OE-OI 
,26ii05806E-01 
.32«712«2E»>01 
.3919<4110E-01 
.«fc58«525E-01 
.5U899a90E-0l 
,63390S59E-Ol 
.7280001/4E-01 
,a2R<»«580E-01 
,9358(hTfe9E-01 
.10a88595E+00 
.116R0199E+00 
,129?9a59E+00 
,ia?.34l26E + 00 
.1559g027E70n 
. 17000928E+00 
.18458558E+00 
.199826S5E+00 
.21510730E+00 
.23l00353e+00 


-.50912247E-01 
-,a 99 Ua 234 E -01 
-.46786989E-01 
-.47342858E-01 
-■. 4 Efe 21 flS 9 E -01 
-.435U276UE-01 
-‘.41078n78E-01 
-'.3B2(411fe5E-Ol 
-.35n983UE-01 
-.S17393B6E-01 
-’.28295190E-01 
-■.24921795E-01 
-. 2174 a 095 e -01 
18825589E-01 
1618flU76E-ni 
-.13830033E-0r 
-‘.1 t724098E«»01 
-,982992fc9E«02 
-’.81 16B513E-02 
-.85672286E-02 
•.5129BB76E-02 
-.37593524E-02 
-.2a200189E-02 
l088l96feE-02 
.256065398-03 
.16331177E-02 
.30512413E-02 
.«5l68636E-02 
.60128731E-02 
.75456025E-02 
.91967362E-02 
.1 S058310E-01 
,i32«3320E-0l 
.15829109E-01 
.186U7723E-01 
'.2U90743e-0l 
.24332496E-01 
,27169B98E-01 
.29979690E-01 
.32745S74E-01 
,35449598E»01 
.38073872E-(U 
.40603093E-01 
.43025607e-0l 
.45336855E-01 
.47539130E-01 
,a9640«95E«01 
,51651U64E-01 
.535e005RE-0l 
'.55a31923E-0l 
.5720961UE-011 
.5fl912a51E-01 
.60540566E-01 
,620926O6E-0i 
655660 14E-01 
.64958730E-01 
.66268100E-01 
.67U92208E-01 
.6B629410E-01 
.69677285E-01 


0BIGI3Sr.AIi BAC4B IS 
OE POOR QUALITY 


6l 


,2a7290l5E+00 

.70633822E-01 

,2639«1UBE+00 

.71496U77E-01 

,2fi09Mii8E + 00 

.72283780E-01 

.29823573E+00 

.72935771E-01 

,3158211 lE+00 

,7350a81UE-01 

.333&662aE+00 

.73975007E-01 

,3517illfc9E + 00 

,74342909E*01 

,370ni929E+00 

.74607337E-01 

,38ftit7050E + 00 

.747869b«E-01 

,<»070B6£ieE + 00 

,74620595E-01 

,U257783iE+00 

.7U766077E-01 

.UUU57722E+00 

.746023688-01 

,ii<>3«3«71E + 00 

.74328132E-01 

.US232163E+00 

,73942521E-01 

.50120B59E+00 

.7344iJl24E-01 

,520fifcfc76E+0Q 

.72830944E-01 

,53B6(S823E + 00 

.72101315E-01 

.55758506E+00 

.71254024E-01 

,5781 8947E+00 

.702875588-01 

.59465457E+00 

.89199941E-01 

.fel29558t>E+00 

.67989380E-ni 

.631 06267E4>00 

.88654380E-01 

.6U8956UOE+00 

.65193S70E-01 

.66661201E+00 

.8360a7nOE-Oi 

.6Ba00«27E4-00 

.818873908-01 

.7031250fiE+00 

,800U1747E-01 

,7l79U2a2E+no 

.58088930E-01 

.734Uii279E600 

.55970O45E-O1 

.75060B17E700 

.55746B0BE-01 

,766a2534E+00 

.51401998E-01 

.781B7389E+00 

.4fl938802E-01 

.79fc94R00E+00 

.48359598E-01 

.811839188+00 

,43fe8798flE-01 

,8259ua55E+00 

.90870923E-01 

,8398648fcE+00 

.37981892E-01 

.B5341131E+00 

,350255498-01 

,88658«57E+00 

.32040020E-01 

.87937578E+00 

.29074484E-01 

.99175958E+00 

.28181554E-01 

.90369481E+OO 

,25408008E-01 

,915t313ftE+00 

.2078798RE-01 

.92801800E+00 

.1B343501E-01 

.93850612E+00 

, 16087738E-01 

.94595079E+00 

. 14O28140E-01 

.95491055E+00 

^121879758-01 

.98314744E+00 

.10507215E-01 

, 970628978+00 

.9045515BE-02 

.97731771E+00 

.77728876E-02 

.983191 42E+00 

,888888398^02 

,9B«22528E+00 

,57854833E-02 

,9923910flE+00 

.50559988E-02 

,99587a49E+00 

.44944813E-03 

,99605372E+00 

,409575108-02 

.99950735E+00 

,3R582107E-02 

.99999989E+00 

9 

,377588B4E-02 

3 2 00,0 

UPPER 8UREACE CALCULATION USING LfC 13,0 PCX THtCK AIRFOIL RESIGNED BY ALLISON 
35, -O. - 

t 2 
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.17762559E-02 
,3tl?5ll 07E-02 
,Sfe829827E-02 
,8563575hE-0? 
.12071072E-01 
.162n803E-0l 
,20<}fl830PE-01 
,2fe«05808E-Ol 
,52U7l2«2E-0t 
.3919ill10E-01 
,U658a5a5E->01 
.5U609««?0E»0l 
,<>3390359E-01 
,728000nE-0l 
,fl2R6«580E-0l 
.93568769E..01 
,] 0Oflfl595E + O0 
.U68fll99E'S-00 
.12929«59E'fOO 
,)a25iJ1£6E + nO 
,15592027E+On 
, 170r>n92RE + 0f) 
.tBU5fi558E+00 
. 1998P635E+00 
,21510730E+no 
,25100353E+no 
.29729015E+00 
.2&39«laRE+00 
,28fi931iiftE+00 
,29eP3373E+00 
.315821 1 1 E+00 
,353fc6629E+00 
,3517ill69E + 00 
,370ni929E^00 
,388i)7050E + 00 
,«0708fcii8E + 00 
,il2577031E+OO 
,ai{/l57722E+0 0 
,«83a397lE+00 
.a6232163E+no 
.50l2ft859E + OO 
,52006676E>00 
.53886823E+00 
,5575850feE4-00 
,576189a7E+00 
,59a65a57E4'00 
.612955S6E+&0 
,fe5106287E+00 
,fc«895fe/iCIE + 00 
.88fe81201EtOO 
.68900B27E400 
.70112508E+00 
,717942fl2E+00 
,759«4279E+00 
,75Ofc0817E+Q0 
.7fc<.«235UE>00 


087« 

.23873T09E+00 
.3879078BE+00 
,520323fe9E+00 
,fe2P0fe53«E+00 
,7?.fi61R3«E + 00 
.79895«59EfOO 
.88S97527E+00 
.9259R157E+00 
,977100fl8E-t-00 
.10228199E+01 
.1 OB97190E+01 
^ in8R2279E+0l 
. 1 1 06577«E+ni 
.11 1 71721E-^01 
,il215306E4'01 
,1 1220805E+01 
, 1 1 19B619E + 01 
,11172711E+01 
.111 29790E+01 
.110880UBE+01 
, llCit2«33E + 0l 
, 1O999997E+01 
. 109556fl4E+01 
, 1 09l53n9E+01 
, loa71666E«01 
.1 9P3a096E<-ni 
, 1079U909E+01 
.1076n5flfcE+01 
. lf>72a985E + 01 
.10693639E401 
. 10P60979E+01 
.in631555E+01 
. 10B999U1E^01 
,1057niBlE+01 
. 10537295E+01 
, lfl506556E + 0l 
. 1007«6?2E+01 
, 10U97203E+01 
. 10U19195E+01 
aO39371feE701 
.1038T261E+01 
.103435362+01 
,103153S2E+01 
.10282740E+01 
.10243370E+01 
. 10202960E+01 
, 10156384E+01 
.10104312E+01 
.lf)OOU108E+Ol 
.99778119E+00 
.99O03184E+00 
.9B1314«BE+00 
.97149016E+00 
,96091 147E+00 
.949196442+00 
.93645976E+00 
.92251687E+00 



,7BlftTS«9E+00 
,79^>9««flOE + 00 

,flP59Ua55E*00 
,fl59Rfc«RBE+00 
,R5%un51E + 0 0 
,R<)R:^R<i57E+00 
.B79J757BF. + 00 
,«917595»E+00 
.903fc9«8iEf00 
.91513138E400 
.92hOieOOE+00 
.93t>3Q«>12E^OO 
,«>a595n79E+on 
.95iJ9l055E + 00 
.9R51«7«UC>'^0 
,970fc2<>97E^00 
.977517718+00 
.9R319142E+00 
.98R2232hE+00 
,992591 08E+00 
,995R7ilii8E + 00 
.99en557?E+00 
.999507356+00 
.999Q99R9E+00 
,1 -.00001? 


,9070Ut'95E + 00 
.B907175PE+00 
.B72085fcRE+00 
.85100503E+00 
.827<>«i929E + 00 
.802?i|672E + 00 
.7770mfe7E + 0fl 
.752S8171E+00 
.75n53U9aE+00 
.71 139979E+00 
.69?03«62E+00 
.680R7085E+00 
.6Rfia8«5tE+00 
,8?7a61?UE+00 
.fcil7591fcfeE + 00 
.fcS855U72E+00 
.830210276+00 
.82231996E+00 
.falii809fe8E+00 
.80678732E+00 
.5983928RE+00 
,5887fl023E+00 
.S78R1908E+00 
.558U7992E+00 
.55518708E+00 
«9 



I(-PUT + SO^^E CnfPUTED OlUlvTITIES 


*•*.*•*. *.*»*p*-*»*»*.*.»ii*-*-*«*.*-*-*»*»*-*-*-*-*-*-*-*-*-*-*-*»*-*»*^*-*-»»*»*-*-*-*-*-*-*-***"*-*«*-*-*-***-»"***"*-*"»"*"*"** 

ITS* so J« 100 TOL* .00001 0Z» ,05000 DS« ,01000 USTgP* ,05000 --F* 1,00000 


IBLCSl 

DISCONTINUOUS SUCTION OR INJECTION RIVE-i BELP“ 




SUCTIC1N OR injection VELOCITY S LOCATION MSI AX* 7 




WWALUM8) 


SOS(HS) 





-3,00000 


,40000 





-5,50000 


.50000 





-a.oonno 


.60000 





-a, 50000 


.70000 . 





-5,00000 


,80000 





-5,50000 


.90000 





-6,00000 

10 

,00000 




INCO 0 

AH1NP3P* ,88506200E400 GAMMA* , 1 UOOOOOOE+O 1 



COMPRESSIBLE ELOWCINCP05, STEWaRTSON 

TRANSEORmaTION used 



B« suction# transition 

ANaVLSIS QE yawed wing laminar 

BPIINDARY LAYER 



INT5« 

1 CM* 1,00000 





lev* 

0 INT9* 161 

RhO* ,3000000 

OE-05 




IBP* 2 







N 

XA(N) 

ZA(N) 

STHCN) 

THIN) 

FSTHIN) 

EZ7H(N5 

1 

-, 10000000E401 

o'. 

0 , 

o'. 

,501117792400 

.32173472E-01 

2 

-,*99«985iet00 

■.318475865-04 

.50250024E-03 

,447916892-01 

.500534862400 

.308960182-01 

S 

•.99B01039ETOO 

.12UT2008E-05 

.19934782E-02 

.892396982-01 

,498856772400 

,277692522-01 

0 

-,9955556aEt00 

,26<><J63a7E-03 

.04525«!69E-02 

.133431032400 

.496177562400 

,23755n57E-01 

5 

-.R921S990E400 

.tt6t624S4E-03 

,787U1690E-02 

'.177555182400 

.492522422400 

,190158952-01 

6 

-.9R776389Et00 

j69008430E-03 

,122556892-01 

.221688032400 

,487895782400 

.133809162-01 

7 

-.98292976E400 

,944515672-03 

,175958342-0! 

.265888102400 

.482288942400 

.625485022-02 

8 

-.97613759E+00 

.12111269E-02 

.238936502-01 

',310191352400 

.475716852400 

-.292533772-02 

«■ 

-,968889122400 

.14722551E-02 

,311468222-01 

,354620762400 

.46826032E400 

-.155540412-01 

10 

•,9606«9S7E400 

.170708012-02 

.39349934E-01 

,399183262400 

.459991442400 

-.248665292-01 

11 

-,05159673E400 

.18927079E-02 

,484944582-0! 

.443877092400 

.45095145E400 

-,370560372-01 

12 

-,9Uia7283E40P 

.a00Ul255E-0? 

.58568974E-01 

.4BB69631E400 

,441194492400 

-, 504730532-01 

13 

-,930«8302E400 

,?01365P5E-02 

.69558788E-01 

‘.533630522400 

,430809232400 

-.650551562-01 

10 

-,91859B03E400 

. 18915196E-02 

,B144a4Q6E-ni 

’.578660352400 

.419851862400 

-.793365032-01 

15 

-.90589329E400 

.1607439UE-02 

,942025596-01 

'.625762632400 

,408406962400 

-.950465822-01 

16 

-.892206401400 

, 11289977E-02 

.107805622400 

.668912062400 

.396611322400 

-.111956862400 

17 

-.877850502400 

.42313630E-O3 

.122220B1E400 

.714072502400 

.384343902400 

-.125914842400 

18 

•,B62689t5E*00 

-,538795U3E-05 

.137ai2fcUE40n 

.759212522400 

.371600532400 

-,139186862400 

19 

-.84679355E400 

-■.17910643E-02 

,153357692400 

■.R0434018E400 

,360330202400 

-.173400042400 

20 

-.83019369E400 

-.34066194E-02 

.170035762400 
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,20ll297RE+Ol 

.59H0R22nE+ni 

-.aR91P6M9£+0n 

,106963708+00 

155 

.9431Rl<iaEtOO 

,66R88839F"02 

,201T2707E+Ot 

.602315758+01 

-,a9388Tu7E+uO 

.102950598+00 

156 

.98B22326E+00 

,57S5ap,33E<‘0? 

,20??3830E+f)T 

,6Ofe57156E+0' 

»,a97P7P3lE+00 

.978659958-01 

1ST 

.99239108E+00 

.50SS99P6E-02 

.2n266101E+or 

.610P50S1E+01 

».50iiao70E+00 

.932135108-01 

158 

.99567«a9F*00 

, aaoaafci 3 E«q? 

,20299U52E+0l 

.615155298+01 

-.503665098+00 

,l!903n6«E-01 

159 

.99Sn5372E+no 

. 009575 loe-os 

.20323576E+01 

.61909233E+01 

-.50536356E+00 

,b5a867iaE-ol 

160 

.99950735E+00 

.38562107E-C? 

.203383098+01 

,623B790?E+01 

-,5065059aE+00 

,82a89585E-01 

161 

,999999ft9E+00 

.37753ft8UE-oa 

,20Sa3?99E+01 

■.62R31B53E + P1 

-.50636002E+00 

.61016105E-01 

IFPT* 

3 1NTI« 2 

NLISTa 0 DXa 0,00000 




CB UPPER 

SURFACE calculation USING LFC 

13,0 PC7 THICK 

AIRFOIL DESIGNED 

BY ALLISUN 


P8I» 35 

,00000 OTRIP* 

0,00000 





1N7RL« 1 

IFR- 2 







RNLCIJ , 

11000000,0 

INTV* 1 L» 83 

fiOU IS COMPRESSIBLE aNO THE FOLLOWING UNCI) 18 THF macm NO. nISTRiHUTION normal TO TmE LtApiNG EDGE (MEN) 


I 

XV 1 1) 

UM(I) 

2 

0,00000 

,08740 

3 

,00065 

.23874 

4 

,00178 

.30791 

5 

.00192 

,52032 

6 

,00568 

,62806 

7 

,00856 

,72062 

8 

,01207 

,79895 

9 

.01621 

.66598 

10 

,02099 

,92548 

11 

.02641 

,97710 

12 

,032a7 

1.02282 

13 

.03919 

1,05972 

la 

.04658 

1.08823 

15 

,05465 

1,10658 

16 

.06339 

1,11717 

17 

,07280 

1,12131 

18 

,08286 

1,12208 

19 

.09357 

1,11986 

20 

,10489 

1.11727 

21 

.11680 

1,11298 

22 

,12929 

1,10884 

23 

.14230 

1.10424 

2 a 

,15592 

1,10000 

25 

.17001 

1.09557 

26 

,18459 

1,09153 

27 

,19963 

1,08717 

28 

,21511 

1.08341 

29 

.23100 

1.07949 

30 

,24729 

1,07605 

31 

.26394 

1.07250 

32 

.28093 

1 ,06936 

33 

.29823 

1,06605 

3a 

,31582 

1,06316 

35 

.33367 

1,05999 

36 

.35174 

1.05702 

37 

.57002 

1 ,05375 

38 

,58847 

1,05066 

39 

.40707 

1,04746 

ao 

,42578 

1.04472 

41 

,44458 

1 ,04192 


\ 






42 

.46543 

1,03937 

43 

.48252 

1,03673 

44 

.50121 

1.03435 

45 

.52007 

1,05154 

46 

.53987 

1.02827 

4? 

.55759 

1.02434 

48 

,57619 

1.02030 

49 

.59465 

1,01564 

50 

,61295 

1.01043 

51 

.63106 

I , 0 0 4 # 1 

52 

.64896 

.99778 

53 

,66661 

.99003 

54 

,68401 

.98151 

55 

.70113 

.97149 

56 

,71794 

,96091 

57 

,73444 

,94920 

58 

.75061 

,93646 

59 

,76642 

,92252 

60 

.78187 

,90745 

61 

, 79695 

,89072 

62 

,81164 

.87208 

63 

,82594- 

,85101 

64 

,83986 

,82767 

65 

,85541 

.80255 

66 

,86658 

,77704 

67 

,87938 

,75258 

68 

,89176 

,73053 

69 

,90369 

,71140 

7D 

,91513 

,69503 

71 

.92602 

,68067 

72 

,93631 

,66848 

73 

,94595 

,65746 

74 

,95491 

.64759 

75 

.96315 

,63855 

76 

.97065 

.63021 

77 

.97752 

.62232 

78 

,98519 

.61465 

79 

,98822 

,60679 

80 

,99259 

,59839 

81 

,99567 

,58878 

82 

,99805 

,57682 

S3 

,99951 

.55846 

84 

1,00000 

.55517 


MGR*D« ,96lOf)OOOE + 02 

VGRAOs 8«,8fct.77 XV(U« 15a<)OOQOE-Oa 


NON-DJMCN8IOMAL DISTANCE EROM LOWER 8UREACE TRAILING EDGE (IF UPPER SURFACE IS TO SE COMPUTED, 18P«i) , OR FROM UPPER SURFACE 
trailing EOSE (IF LOWER SURFACE IS TO RE COMPUTED, I SPrO ) TO ATTACHMENT LINE ■ 8ATT/CH ■ , 1 OlZaSMSE+Ol 


SUCTION, transition ANAYLSIS OF YAWED WINC LAMINAR BnuNOARV LAYER 


UPPER SURFACE CALCULATION USING LFC 13.0 PCT THICK AIRFOIL DS8ICNEO PY ALLISON 


REYNOLDS NUMBER DEFINED PY RHb OL/NU 

VELOCITV ORADIEFT AT ATTACHMENT LINE* S8.67 


XV 

,lSA90fl00Ea01l 

|fc 50302 a 7 E -03 
,177fc2519e-02 
.31I251107C-02 
,56«29«27C-02 
.65fc357?feE-02 
,1207l072e-0j 
ISnSE-Ol 
,2O9803OOE»O1 
.26il05l»06E-01 
,32fl71?42E**01 
,391901 tOE-Ol 
.a<)5B«525e-01 
,5Ofc«9O90E»Oi: 
,63S90359E-01 
.728000i4E»01 
,82864580E»ni 
.9356fc769E-tVt 
,100B8595€+00 
. llfePf)}99E*00 
,12929U59E+t)0 
.ia 234 iaiEtOP 
.i5592027E>00 
, 17nn0928E+00 
, 184585?8E*00 
.199828396+00 
,21510730E+0(V 
, 231033552+00 
.2U729015E+00 
,263941482+00 
,28093146E+00 
,298233732+00 
.S1582111E+00 
.33366824E+00 
,351741692+00 
,370019292+00 
,188470502+00 
,407066482+00 
,425778512+00 
,444577222+00 
,0634+4712+00 
,482321632+00 
,501208592+00 
,520066765+00 
,538868232+00 
,557585062+00 
,576189472+00 
,594654372+00 
,612953862+00 
,631062672+00 
,648956402+00 
,666612012+00 
,684008272+00 
,70112508E+-eb 
,717942822+00 
.75444279E+00 
.750608172+00 
,766423342+00 
,781873892+00 
.796948002+00 
,8l 1638182+00 
,825942552+00 
,839864862+00 


SXV 

0 . 

,15752780E-P2 
,42«15860E+02 
,727&77162»02 
,105623312-01 
,141931 13E-01 
,182364602-01 
,22730621E-0V 
.277101552-01 
,331989322-01 
.392185822-01 
,457902242-01 
.52936238E-0J 
,606796292-01 
,690398732-01 
,780297652-01 
,876519272-01 
,97R99»t'082«01 
,lQ87s"832+00 
,120218762+00 
.132255852+00 
,144654102+00 
,157992762+00 
,171651482+00 
,185809152+00 
,200444152+00 
.215534652+00 
,231057312+00 
,246988n4E+00 
,263302712+00 
,279976388+00 
,296983702+00 
.514298922+00 
,331895572+00 
,349746892+00 
,367826082+00 
,386105602+00 
.40455750E+00 
,425153552+00 
.441865062+00 
,460665082+00 
.479524572+00 
,498415422+00 
.517308962+00 
,536177102+00 
.554992722+00 
,573726722+00 
,592358212+00 
.610855122*00 
,629194602+00 
,64735256E+fl0 
,665305832+00 
.683032772+00 
,70Ci51359E + 00 
.71772962E+0Q 
.734662672+00 
,751295602+00 
.767613152+00 
,783601202+00 
.799246892+00 
,814540072+00 
,629474762+00 
,844050032+00 
,858268942+00 


SXVINC 

0 . 

,254209702-02 
,621925372-02 
.1Q592972E-01 
,145536192-01 
, 187049732-01 
.236015262-01 
,273316712-01 
.317516972-01 
.3626706UE-01 
,408837992-01 
,456122462-01 
,504745262-01 
,555100112-01 
,607701592-01 
,663111902-01 
,721796182-01 
,784075982-01 
,850154152-01 
,920103102-01 
,995950502-01 
,1 07170672+00 
, 1 15330432+00 
.123866612+00 
,132769192+00 
,142027152+00 
.151630002+00 
.161564442+00 
.171815122+00 
.182367102+00 
.191204102+00 
.204309642+00 
.215667122+00 
,227259102+00 
,239068342+00 
,251079282+00 
,263276002+00 
,275641872+00 
,288158242+00 
.300804072+00 
,313557572+00 
,326398632+00 
.339306812+00 
.352261122+00 
,365243862+00 
,378244152+00 
,391252842+60 
. 404257682+00 
,417246062+00 
,430209762+00 
.443142252+00 
.456037542+00 
,468892742+00 
,4817P»43E+60 
.494484992+00 
,507220292+00 
,519912142+00 
,532560742+00 
.545166462+00 
,557730002+00 
.570256242+00 
.582757732+00 
.595255182+00 
.607773192+00 


II 

0 . 

,987501442-01 

,269740552+00 

,438283282+00 

,587895132+00 

,709626382+00 

, 814200972+00 

,902710312+00 

,978434582+00 

,1 04566862 + 01 

,110399142+01 

,115564722+01 

.119733872+01 

, 122954932+01 

.125028242+01 

.126225292+01 

.126695152+01 

.126779882+01 

, 1265?9202*01 

. 126236482+01 

.125751532+01 

.125284422+01 

.124764512+01 

,124285042+01 

,125784372+01 

.123328162+01 

.122835082+01 

, 122410592+01 

.121967832+01 

, 121579122+01 

.121177782+01 

.120823612+01 

,120448952+01 

,120122152+01 

,119764952+01 

.119426572+01 

.119057142+01 

. 118709852+01 

,118349022+01 

. 118039222+01 

.117722772+01 

,117434892+01 

,1 17135992+01 

,116667932+01 

, 1 16549482+01 

.116181012+01 

,115736192+01 

.115279612+01 

,114753362+01 

,114165022+01 

. 1 13484802+01 

.112735742+01 

,111860172+01 

.110675222+01 

.109765212+01 

.106569962+01 

.107246322+01 

,105807252+01. 

.104231892+01 ' 

.102529202+01 

,100639002+01 

.905316262+00 

.961520242+00 

,935154022+00 


THXV 

,3133721?E+01 
,314159262*01 
.319260022*01 
,322590892+01 
.325«b742f. + 0l 
.329250692+01 
,'32693702+01 
.336177402+01 
,339693692+01 
.343236272+01 
.346803752+01 
,350396772+01 
.354017592+01 
,357668542+01 
,361350402+01 
,365061932+01 
.368799402+01 
,372558002+01 
,376333222+01 
,380120782+01 
.383917722+01 
.387721862+01 
.391531072+01 
.395343972+01 
,399159492+01 
.402976802+01 
,406795342+01 
,410614432+01 
,414433472+01 
.418252072+01 
,422069822+01 
,425686382+01 
,429701462+01 
.433514702+01 
,437325822+01 
.441134582+01 
,444940672+01 
.448743762+01 
.452543482+01 
,456339492+01 
,460131522+01 
,463919312+01 
.467702492+01 
.U7'(4S062E+01 
,475255332+01 
,479020392+01 
,482781522+01 
,486536392+01 
.490284712+01 
.494026372+01 
,497761292+01 
,501489442+01 
,505210992+01 
. ,508926392+01 
%^5i 2636232+0 I 
.Ir5l634iu2+01 
-:i««2004 1762 + 01 
.‘•-S23739 282+01 
"^527435002+01 
-■‘',531 130532 + 01 
.534826172 + 01 
,538531442+01 
.542245672+01 
.545978302+01 


FUTH 

-.554200162+03 
-.534286472+03 
.237624122+03 
-.772115872+02 
-.205949672+02 
-.158347462*02 
146886482+02 
-, 114749792+02 
-.58201736F+01 
-.871300132+01 
-.372985382+01 
-.940932392+01 
-.623647502+01 
-, 100636722+02 
-.573770982+01 
-,605065562+01 
-.177871772+01 
-.335819572+01 
,103850032+01 
-.252458972+01 
,108007882+01 
-,101856852+01 
,635339092+00 
-.632397752+00 
,833473032+00 
-.859288992+00 
.10900608F+01 
-.675171012+00 
,858833212+00 
-.538244852+00 
,770638842+00 
-.606702352+00 
.805741902+00 
-.649262082+00 
.528527242+00 
-.621844622+00 
,514776032+00 
-,456353072+00 
,7360471 12+00 
-,378208502+00 
,485515102+00 
-.364331042+00 
,575144082+00 
-.637514882+00 
-,164053572+00 
-.840303282+00 
.264931432+00 
-.750508072+00 
-.269855292+00 
-.869230612+00 
-.247224172+00 
-,116393462+01 
-.650733592+00 
-.112592692+01 
-.378102352+00 
-,114970062+01 
-,704816922+00 
-.114398152+01 
-.735531532+00 
-.151176982+01 
-.140534042+01 
-.316170142+01 
-.170601032+01 
-.150759372+01 


F3V8INC 
,488721622+00 
.426551502+01 
,968223702+01 
,153933672+02 
,186913752+02 
, 198538542+02 
.207079752+02 
.202451742+02 
.199726232+02 
,193869952+02 
,185246312+02 
.168903992+02 
.139176802+02 
, 1009flOPUE+02 
.599390062+01 
,285604082+01 
,840404892+00 
-.386536652+00 
-.863531842+00 
-,112480562+01 
-.153837942+01 
-.125289072+01 
=,121833682+01 
-.1 12350672+01 
-,103835242+01 
-.996407732+00 
-.919688922+00 
-.825901702+00 
-,770477652+00 
“,702430442+00 
-,650648112+00 
-.611413822+00 
-,570690932+00 
-.538236162+00 
-,536899062+00 
536992902+00 
», 536967242+00 
-.530074792+00 
-,461680902+00 
-,437709332+00 
», 420915072+00 
-.405666272+00 
-.362756572+00 
-,365396972+00 
-,458916952+00 
-.549684302+00 
-, 598351872+00 
-,637830342+00 
-.725491922+00 
-.821225582+00 
-.91198931E+00 
«, 102546522+01 
-,116693382+01 
-.139788852+01 
-.140036532+01 
-.149885212+01 
-.161604312+01 
-.173365112+01 
-.183416762+01 
-.194217062+01 
-.309411062+01 
-.237120352+01 
-,244099122+01 
-.25400752F+01 


,85Sai )3i£4«0 ,872l3a23E+00 
,8fefcSBU57E400 ,8B58ai56E*60 
.87'»37576E400 ,89877203E + 00 
,a<»17!956E*00 ,91ia892fcE + 00 
,90369a8lE*no ,9237a2SlE+00 
.9151313RE+00 .935875382+00 
,9260ieoo€+00 ,<?866330aE + 00 
.93fe30612E+00 ,957lfc55fcE+00 
,9a595079E+00 ,987027fe9E+0D 
,95a9l055E+00 ,9781T85iE+00 
,963ia7aaE400 ,9Ba56HfeE+00 
,97062h97E+t)0 . 99220256E + 00 
,97731771t-+00 ,99901292E + O0 
,985l9ia2E+00 , 1 00U9B5flE + 0 1 
,9S822328E+D0 . 1 0 1 0098 lE+0 1 
,992391082+00 ,10ia3293E+01 
,99567aaBE+00 .101766032+01 
.99805372E+00 , 1 0200727E+0 1 
.999507S5E+00 . 10218860E+01 
.99999989E+00 ,102200495+01 


.62032907E+00 
.63291601E+00 
,6a5a9958E+00 
.65799569E+00 
.67030076E+00 
.6B23015aE+00 
,69389219E+Q0 
,70a97925E+00 
,715080672+00 
.725323772+00 
,730000172+00 
.702780812+00 
.750290602+00 
.7S692BttlE+00 
,762606892+00 
,76701«672+00 
,77) 198502+00 
,77396093E+oO 
,775676582+00 
.776261602+00 


. 90+768932+00 
.877951682+00 
,850315202+00 
.B25UO53S2+O0 
,803785262+00 
. 7652914802+00 
.769291702+00 
.75529682E+00 
.702802352+00 
.731690732+00 
.721080212+00 
.7120521 22+00 
.703137102+00 
,69U07076E*00 
.685587362+00 
.676)02562+00 
.665201792+00 
.651727312+00 
.651006552+00 
.627263092+00 


,509738502+01 

.553535952+01 

.557378332+01 

,561269832+01 

.565210302+01 

.569196792+01 

.573225101+01 

.577291002+01 

.581392)82+01 

.585520172+01 

,589680672+01 

.593871352+01 

,598082202+01 

.602315752+01 

.606571562+01 

,610850512+01 

,615155292+01 

,619092332+01 

,623879022+01 

.626252202+01 


-.362706032-01 
,) 00696172+01 
,221715912+01 
,239571332*01 
, 221100822+01 
,160695072+01 
,)3325073E+0i 
,937050322+00 
.863662772+00 
,500222372+00 
,517079312+00 
,297893282+00 
,217906202+00 
-.172080532+00 
-,966592612-01 
-.125505192+01 
.833238002+00 
-, 102796902+02 
,180753752+02 
-.903768762+01 


,251665132+01 
,230701002+01 
.207061052+01 
,176038602+01 
,107872802+01 
,126253702+01 
, 1 10818002+01 
,999789782+00 
.923953212+00 
.877232002+00 
,658716272+00 
,865076992+00 
,908719102+00 
,102203002+01 
,121030332+01 
,166951822+01 
,219168162+01 
,537325322+01 
.801310822+01 
,206210222+01 
















0U7PU7 


© 


**+ LEiOING- 

EDGE CONTiMlNtTION 1287 





REYNOLDS NUMBER* IIOOOOOO, RTHETA* 51,2 

)(■ -.000015 8* 0,000000 SCOMP* 0,000000 

U» 0, 

NO TURBULENT CONTiMlN+TION AT A.L, 

000000 AHE3D« ,082902 0U/D(S/L3* 72.795397 12 ITERATIONS 


DELTAl* 

0750 TH2T61* .1975 

COU/DZ5Z90 

■2.2006 

(DV/DZ)Z«0 ■ .1102+01 airfoil SLOPE" .900E+02 


Z 

ZCOMP 

w 

U 

V 

8TV CFV T 

RHOO 


0,000000 

0,000000 

-.289171 

O.OOOOQO 0 

.000000 

1.156667 

1,226681 


.050000 

.058772 

-.290886 

.068570 

.038029 

1,156592 

1,226760 


.100000 

.077539 

-.295902 

,133679 

.075508 

1,156373 

1,226993 


,200000 

.155000 

-.515399 

.253793 

.108809 

1.155526 

1,227893 


.300000 

.232067 

-.306216 

.360975 

.219829 

1.150176 

1,229328 


,000000 

,309788 

-.387137 

,055969 

.288029 

1.152379 

1,231206 


.500000 

.386975 

-.036980 

.539592 

.350020 

1.150195 
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Figure 3. Upper surface skin-friction distritutions for sample case 


